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Zur neueren Entwicklung der Schalentheorie*) 
Von W. ZERNA 


1. Einleitung 


In den letzten Jahren haben Schalentragwerke eine umfassende Bearbeitung erfahren. Die 
Anzahl der Veröffentlichungen in aller Welt dürfte die Zahl 1000 bereits weit überschritten haben. 
Angeregt wurde diese intensive Forschungstätigkeit durch die vielen Anwendungen, die Schalen 
auf den verschiedensten Gebieten des Bauingenieurwesens, des Maschinenbaus, des Flugzeug- 
baues, ja selbst der Raketentechnik und Kerntechnik haben. Schalen werden heute aus allen 
praktisch vorkommenden Werkstoffen hergestellt, aus Stahl, Beton, Leichtmetall, Kunststoffen 
und dergleichen. 

Die Schalentheorie läßt sich in drei Teile aufgliedern: Statik der Schalen, Dynamik der 
Schalen und Stabilitätstheorie der Schalen. Jedes dieser Gebiete stellt wiederum einen umfas- 
senden undin viele Einzelfragen verästelten Wissenszweig dar. In den folgenden Ausführungen 
soll nur die Statik der Schalen behandelt und insbesondere auf die Entwicklung der allgemeinen 
Theorie eingegangen werden, da hier in den letzten Jahren besonders beachtliche Fortschritte 
erzielt worden sind. 

Für die Entwicklung der allgemeinen Theorie ist die Einführung des Tensorkalküls in der 
von Rıccı angegebenen Schreibweise von ausschlaggebender Bedeutung geworden. Damit wurde 
die Möglichkeit geschaffen, eine allgemeine Schalentheorie aufzustellen und die äußerst verwickel- 
ten Zusammenhänge in übersichtlicher Weise zu schreiben. Es sind hier in den Jahren von 1940 
bis 1950 eine Reihe von Arbeiten entstanden, die die Zweckmäßigkeit dieses Werkzeuges unter Be- 
weis stellen. 

Ausgehend von einer solchen allgemeinen Theorie lassen sich dann die beherrschenden 
Gleichungen für irgendeine spezielle Schalenform und für jedes gewünschte Koordinatensystem 
als Sonderfall auf rein algebraischem Wege finden. Darüber hinaus lassen sich dann aber auch 
mit Hilfe allgemeiner übersichtlicher Theorien neue Einblicke in grundlegende Zusammenhänge 

ewinnen. 
: Neben der klassischen linearen Theorie mit ihrer Beschränkung auf infinitesimale Verfor- 
mungen und einem linearen Spannungs-Dehnungsgesetz sind auch Anfänge für nichtlineare 
Theorien zu finden, die entweder geometrisch oder physikalisch oder in jeder Hinsicht nichtlinear 
sind. Die folgenden Ausführungen beziehen sich jedoch lediglich auf die im Sinne der klassischen 
Elastizitätstheorie lineare Schalentheorie. 


2. Gleiehgewiehtsbedingungen und Spannungsfunktionen 

Der allgemeine Aufbau der Schalentheorie erfolgt derart, daß durch Integration der Span- 
nungen über die Schalendicke Schnittgrößen als Funktion zweier Flächenkoordinaten eingeführt 
werden. Es ergeben sich der Längskrafttensor n°%, Momententensor m*?, Querkrafttensor q* 
(griechische Indizes durchlaufen die Zahlen 1, 2). Die Forderung, daß alle an einem Schalen- 
element angreifenden Schnittgrößen und der Belastungsvektor p* im Gleichgewicht sein sollen, 
liefert ein System gekoppelter partieller Differentialgleichungen für die Schnittgrößen, das in 
Tensorschreibweise folgende Gestalt hat: 


bet; 
nPbuatQatp=0; | 

mb, — gP —(, ( 
Eu (m? — b, meP) = 0 | 


*) Hauptvortrag, der auf der wissenschaftlichen Jahrestagung 1960 der GAMM (19.—23. April 1960 in 
Freiberg/Sa.) auf Einladung der Tagungsleitung gehalten wurde. 
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Darin bedeuten bs die Koeffizienten der zweiten Fundamentalform der Schalenmittel- 
fläche. Vertikale Striche deuten kovariante Differentiationen an, und &ug Ist der in der Tensor- 
rechnung bekannte schiefsymmetrische Tensor. Die letzte Gleichung in (1) stellt eine Identität 
dar, die lediglich die Symmetrie des Spannungstensors zum Inhalt hat, so daß also insgesamt 
fünf Gleichungen des Gleichgewichts vorhanden sind. Durch Elimination der Querkräfte lassen 


sie sich auf drei Gleichungen reduzieren. Diese lauten: 
ee: RE, (2). 
n°P bug + ml +p = _ 
Die Gleichungen (2) sind die notwendigen und hinreichenden Bedingungen für das Gleich- 


ewicht. 
: Wird der Belastungsvektor Null gesetzt, so lassen sich die Gleichgewichtsbedingungen (2) 
identisch durch den Ansatz 
ne — er A ++ RD, N ea... 
m — er — du Pre J 


erfüllen. Darin stellen %,, Y, ® vier Spannungsfunktionen dar, die sich als Komponenten zweier 
Vektoren 
ee a 


Y-WUHVN, 


deuten lassen (X, = Grundvektoren der Mittelfläche, Y, = Normalvektor). Zwischen ®* und Y,, 
. YW besteht die Beziehung 

| DP = EP* Yoa+b&YP); 
so daß vier unabhängige Spannungsfunktionen übrigbleiben. Die in (4) eingeführten Vektoren 
lassen sich anschaulich als resultierende Kräfte und Momente deuten, die längs eines Bogen- 
stückes auf die Schalenmittelfläche von der einen auf die andere Seite übertragen werden. 


3. Verformungszustand 


Der Verformungszustand der Schale wird durch den dreidimensionalen Verzerrungstensor 
yır (uk=1,2,3) beschrieben. Der Verzerrungstensor läßt sich durch den dreidimensionalen 
Verschiebungsvektor V,; ausdrücken. Mit Hilfe des Hooxeschen Gesetzes werden dann Beziehun- 
gen zwischen dem Verzerrungstensor und dem Spannungstensor r’* hergestellt und dann schließ- 
lich durch Integration der Spannungen über die Schalendicke der Zusammenhang zwischen 
den zweidimensionalen Schnittgrößen und den Komponenten des Verschiebungsvektors der 
Mittelfläche gefunden. Bei diesem Prozeß treten einige charakteristische Schwierigkeiten auf, 
die nachfolgend umrissen werden sollen. 

Wird stellvertretend für irgendeine Schnittgröße S geschrieben und für irgendeine Ver- 
schiebungsgröße V, so führt der oben angegebene Weg auf einen Ausdruck der Form: 


SEVEN 


Darin stellt G eine geometrische Größe dar, die in Form einer Potenzreihe nach Potenzen 
der Koordinate z normal zur Schalenmittelfläche gegeben ist, also: 


G=eg trzn ri re a ee 
Die Verschiebungsgröße V muß mit Rücksicht auf den gewünschten zweidimensionalen 
Charakter der Theorie ebenfalls als Potenzreihe der Form 

V=V,+2:V,+2W,+:-- ee Tr EEE 


angesetzt werden. Die in (7) auftretenden Koeffizienten V,, V) Va... können als die eigent- 
lichen Unbekannten der Schalentheorie angesehen werden. Esist nun nicht ohne weiteres möglich, 
sie allein auf mathematischem Wege zu gewinnen. Es hat nicht an Versuchen gefehlt, eine mathe- 
matisch strenge Theorie aufzustellen, was allerdings bisher als nicht völlig gelungen angesehen 
werden muß. 


Man geht nun so vor, daß im Hinblick auf die geringe Dicke der Schale die Reihe (7) nach 
dem linearen Glied abgebrochen wird. Demgemäß ergibt sich: 


V=u,+uz en Un Lupe Eee BEE 


ı das eh Gesetz in eine ot ST Form zu bringen, in 
weit angenommen, daß die Spannungen in Normalenrichtung Null gesetzt werden können. 
2 ‚entsprechende Komponente des Spannungstensors nimmt damit. den Wert 


. Mit Hilfe dieser Annahmen gelingt es an v, allein Auge Do a ücken und damit en nt 
R. ‚alle Schnittgrößen n*P, m*#* durch v, auszudrücken. . | BE 
En _ Dieses Vorgehen hat aber eine Reihe von Widersprüchen zur Folge. Aus (9) folgt g* = 0. Be; 
_ Weiterhin sind y3; = 0 und 73; = O miteinander unverträglich. Es muß darauf verzichtet werden, 
1 Ei durch v, auszudrücken, was bei der weiteren mathematischen Behandlung nicht slor]: da in 

den drei Gleichgewichtsbedingungen (2) die Größe g nicht auftritt. 

Zu den aufgezeigten Widersprüchen tritt nun noch als besonders unbefriedigend hinzu, 
daß wegen (8) bei der Bestimmung von (5) keine einwandfreie Aussage möglich ist, wieviele Rei- 
_ henglieder mitzunehmen sind. Werden (6) und (7) gemäß (5) miteinander ausmultipliziert, so 
ergibt sich: 7 


E. | S=-9n +2 +) +2HR + HH t+Hd)+" ::...(M. 


=... Da wegen (8) v, verschwindet, muß der Koeffizient bei z? als fehlerhaft bezeichnet werden. 
_ Verläßlich sind im Hinblick auf die Reihenentwicklung nur das erste und zweite Glied in (11). Da 
nun aber v, mit Hilfe der Normalenhypothese bestimmt worden ist, ergeben sich auch Zweifel an 
- der Genauigkeit des linearen Gliedes der Reihe (11); denn es wäre möglich, daß bei Ersatz der 
- — Normalenhypothese durch den genauen paysicahsehen Zusammenhang dies Glied eine andere 3 
Gestalt annimmt. y 

Je nach der Art und Weise, wie jetzt die Entscheidung bezüglich der Mitnahme, des Fort- 
lassens oder irgendwelcher weiterer Spekulationen in der Theorie fällt, ergeben sich verschiedene 
Schalentheorien. 

Das allgemeinste Elastizitätsgesetz, das eine mit der geschilderten Annahme aufgebaute 
Schalentheorie zu liefern vermag, lautet: BE 


n“®= DHPAa, +{BGe o,.}; 
m®P = BHeeA wo, + {BG a, .} 


ET AAN: 


(12). 


Darin stellen D und B die Dehnungssteifigkeit und die Biegesteifigkeit der Schale dar, die 
sich durch die Schalendicke t, den Elastizitätsmodul und die Querkontraktion ausdrücken lassen. 
Hs»ßei und G*?e% stellen Tensoren dar, die von der geometrischen Form der Schale abhängen, 
&,, und w,, drücken sich durch den Verschiebungsvektor der Schalenmittelfläche aus. 

In (12) müssen die in den geschweiften Klammern befindlichen Ausdrücke entsprechend den 
gemachten Ausführungen als fraglich angesehen werden. 

Es ist viel Arbeit darauf verwendet worden herauszufinden, welchen Einfluß die angeführten 
Annahmen ausüben und wie eine in Bezug auf die Größenordnung der einzelnen Glieder einiger- 
maßen einwandfrei aufgebaute Theorie aussehen müßte. Um diese Frage einwandfrei beant- 

-  worten zu können, wäre es erforderlich, eine Theorie zu entwickeln, die keinen Gebrauch von den 
- _ dargelegten Annahmen macht. Dies scheint jedoch bisher in einwandfreier Weise noch nicht ge- 
lungen zu sein. 


nt IN HE an din „aa ale ni, NEN, 
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4. Verträglichkeitsbedingungen 


Die Tatsache, daß sich die sechs Komponenten des Verzerrungstensors durch drei Kom- 
ponenten des Verschiebungsvektors ausdrücken lassen, zeigt, daß die einzelnen Komponenten 
des Verzerrungstensors nicht unabhängig voneinander sind. In der Schalentheorie drückt sich 


dies durch folgende Verträglichkeitsbedingungen aus: 
GP, — 58 He, = 0, 
GP DaB = I: Klin —— 0 » | EEE ‚erbre ia Let ja a Mer/\e 13)3 


&u5 (GP — bE Het) = 0 
i* 
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Darin ist 1 
GR = re IWi.. + (b2 V)ıa Die; bu dar Vorrt » 


Hr = Zen eV ujgt Vern— due W}; 


(V., W = Verschiebungsvektor). 


Die Verträglichkeitsbedingungen zeigen, daß formal zwischen ihnen und den Gleichgewichts- 
bedingungen eine Analogie besteht. Allerdings muß darauf hingewiesen werden, daß der Tensor 
H°® symmetrisch ist, während im allgemeinen der Tensor m*P# nicht symmetrisch ist. Die Analogie 
wird daher nur dann vollständig, wenn auch der Tensor m°P als symmetrisch vorausgesetzt 
werden kann. In vielen Fällen ist dies möglich, und es bestehen dann zwischen den einzelnen 
Größen folgende Analogien: 

GP on, 


HP —meP, 
v>V; 
v_>2W. 


5. Spezielle Lösungstypen 


Die aufgezeigten Schwierigkeiten in der Schalentheorie sind zwar im Hinblick auf den 
mathematischen Charakter der Theorie außerordentlich störend, praktisch fallen sie jedoch weit 
weniger schwer ins Gewicht, da im Anwendungsfall im allgemeinen Sonderfälle vorliegen, bei 
denen sich grundlegende Vereinfachungen ergeben, die die Theorie in der einen oder anderen 
Weise auch in mathematischer Hinsicht verbessern. 

Im allgemeinen wird die Lösung irgendeiner Schnitt- oder Verformungsgröße von der Form 
sein: 


s=S$ (0,, B,, )) : 


Darin sind 6,, 6, die Flächenkoordinaten der Schalenmittelfläche und A = {/L (t = Schalen- 
dicke, L = Längenabmessung). Die unterschiedliche Größenordnung von S und das Verhalten 
bei der Differentiation bezüglich einer Änderung der Größenordnung ermöglicht eine Klassifi- 
zierung in 

a) Membrantheorie, 
b) Lösungen vom Typ der dehnungslosen Verbindungen, 
c) Randstörungen. 


. Ohne auf diese drei Fälle jetzt näher einzugehen, soll festgestellt werden, daß sich bei ihnen 
die aufgezeigten Mängel der Theorie nicht mehr bemerkbar machen. 


6. Membrantheorie 


Werden die Momente und Querkräfte in den allgemeinen Gleichungen des Gleichgewichts 
vernachlässigt, so erhält man die Membrantheorie. Sie beschreibt einen biegungsfreien Spannungs- 
zustand. Es bleiben jetzt nur drei Gleichungen für den symmetrischen Längskrafttensor n%# 
übrig. Das Problem besteht nunmehr darin, festzustellen, ob bei einer bestimmten, vorgegebenen 
Belastung und den Randbedingungen ein Membranspannungszustand möglich ist. Für Schalen 
mit positiver Gaussscher Krümmung läßt sich diese Frage recht allgemein beantworten. Das 


Problem führt hier auf Differentialgleichungen von elliptischem Typ, die mathematisch weit- 
gehend erforscht sind. 


Bei negativer Gaussscher Krümmung der Schalenmittelfläche liefert die Membrantheorie 
unter Umständen singuläre Lösungen bei bestimmten Belastungen. Das hängt damit zusammen 
daß solche Schalen in gewissen Fällen infinitesimale Verschiebungen erleiden können. Aber selbst 
wenn sich endliche Lösungen ergeben, ist die Frage nach den möglichen Randbelastungen nicht 
ohne weiteres zu beantworten, weil jetzt Differentialgleichungen von hyperbolischem Typ auf- 
treten, bei denen keine Randbedingungen im üblichen Sinne vorgeschrieben werden können. 


7. Randstörungen 


a a Fällen ist es möglich, die Membrantheorie als Partikularlösung der allgemeinen 

en geht man dann so vor, daß die Schale zunächst nach der Mem- 
e andelt wir ie ne R { PL 

& wird und die Randbedingungen erfüllt werden, soweit die Membrantheorie 


En /NRR ” wor “ hy Pr , E u ” 
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E: weiterer Randbedingungen. Für den Fall, daß die Randkurve nirgends mit den Aymptotenlinien 

der Mittelfläche zusammenfällt oder sie nur berührt, was bei positivem Gaussschen Krümmungs- 
maß stets der Fall ist, läßt sich eine verhältnismäßig einfache Näherung zur Lösung dieses Rand- 
störungsproblems angeben. Wenn der Rand auf Flächen mit verschwindender oder negativer 
Gaussscher Krümmung mit den geradlinigen Erzeugenden bzw. Asymptotenlinien zusammen- 
fällt, versagt die Randstörungsrechnung in ihrer einfachsten Form. Die Schnittgrößen und 
Verschiebungen klingen nämlich in diesem Falllangsamer ab alsin dem vorher genannten. Schalen 
mit negativer Gaussscher Krümmung wird man deshalb — insbesondere wenn auch das beachtet 

' wird, was bezüglich der Membrantheorie gesagt worden ist, — in allen nicht ganz einfachen Fällen 
von vornherein nach der allgemeinen Theorie behandeln. 


8. Flache Schalen £ 


r Eine flache Schale ist dadurch gekennzeichnet, daß die Tangentialebene in jedem Punkt der 

Mittelfläche eine geringe Neigung gegen die Grundrißebene aufweist. Für derartige Schalen sind 
besondere Näherungstheorien entwickelt worden. Mit ihrer Hilfe lassen sich flache Kugelkalotten, 
hyperbolische Paraboloide und Translationsschalen vorteilhaft behandeln. 

Die Theorie der flachen Schalen ergibt sich unmittelbar aus der allgemeinen Theorie, wenn 
zusätzlich weitere Vereinfachungen eingeführt werden, die mit der besonderen Geometrie der 
Schalenmittelfläche flacher Schalen zusammenhängen, d. h. vorausgesetzt wird, daß die Geometrie 
auf der Schalenmittelfläche euklidisch ist. Damit gilt diese Theorie dann auch für nicht mehr 
notwendigerweise flache Schalen, bei denen das GAausssche Krümmungsmaß Null ist, da diese 
Flächen auf der Ebene abwickelbar sind und daher auch euklidische Maßbestimmungen haben. 

In der Theorie der flachen Schalen werden die Schnittkräfte durch eine Spannungsfunktion 
dargestellt. Die Momente werden durch die Ableitungen der Normalverschiebung w ausgedrückt. 
Diese beiden Funktionen lassen sich aus zwei gekoppelten partiellen Differentialgleichungen 

4. Ordnung bestimmen. Diese lauten: 


7729 —-EiV!w=0, 
21-M)Rg+ELPPw=12Z(1—). 


Hierin bedeutet Z die Lastkomponente senkrecht zur Schalenmittelfläche. V? und V% 
bezeichnen Differentialoperatoren 2. Ordnung. E ist der Elastizitätsmodul, » die Querkontrak- 
tionszahl und f die Schalendicke. Die Lösung dieser Gleichungen läßt sich in vielen Fällen der 
praktischen Anwendung durch Reihen darstellen, deren Berechnung mit Hilfe von elektronischen 
Rechenanlagen mühelos möglich ist. 
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Numerisch stabile Interpolationsformeln 
mit günstiger Fehlerfortpflanzung 
für Differentialgleichungen erster und zweiter Ordnung 
Von ERWIN FEHLBERG 


Für Differentialgleichungen erster und zweiter Ordnung wird ein die Fehlerfortpflanzung charakterisieren- 
der Parameter angegeben. Es werden numerisch stabile Interpolationsformeln günstigster Fehlerfortpflan- 
zung für eine verschiedene Anzahl von rückwärtigen Stützstellen mitgeteilt. 


A parameter is obtained for the error propagation of inter polation formulas for Ist and 2nd order differen- 
tial eguations. For different numbers of backward points numerically stable interpolation formulas of smallest 


error porpagation are presented. 
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dies gestattet. Die allgemeine Theorie dient dann nur noch als homogene Lösung zur Erfüllung | 
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1. Einführung | 
? } za \ 2 ’ j ie 
Zur numerischen Integration von gewöhnlichen Differentialgleichungen sind zahlreic 
Interpolationsformeln vorgeschlagen worden. In dieser Arbeit sollen für die Differentialgleichun- 
gen erster und zweiter Ordnung die möglichen Interpolationsformeln bezüglich numerischer Stabi- 
lität und Fehlerfortpflanzung untersucht werden. Die dazu erforderlichen Kriterien werden, 
soweit sie nicht schon aus Arbeiten anderer Autoren bekannt sind, hier hergeleitet werden. 


Wir schreiben unsere Interpolationsformeln in der folgenden Form: 
für Differentialgleichungen erster Ordnung: 


Mm m R ® 
Ynrı = DAruye—u th & Ar Yu ee 2 2 
= 


yu=— 


für Differentialgleichnngen zweiter Ordnung: 


1 ki A 2 ! 2 [73 [2 
Yr+ı Ze 22 Ar—u Ge —u ZAr—u Y—ut h- = Ar—ı Ye—u» 
N ne Re sa La 
Mm m m , n a 
Yyı= 3 Bu th- 2 Bi Year BA Br—ı Ye—u 
u=0 A= 


Die Einführung der Werte y,_„ auch in die erste Formel (2) unterscheidet unsere Formel (2) 
von den im allgemeinen benutzten Interpolationsformeln für Differentialgleichungen zweiter 
Ordnung. 


Teil 1: Formeln für Differentialgleichungen erster Ordnung 
2. Numerische Stabilität der Interpolationsformeln 


Wir können uns hier auf die grundlegende Arbeit von H. RuTISHAUSER!) beziehen, nach 
welcher die Interpolationsformel (1) stabil ist, wenn das Polynom: 


5 Al. mr —0 (Ay eins Dh a 


u=—1 
außer der trivialen einfachen Nullstelle A = + 1 nur Nullstellen A mit [A| < 1 besitzt. 


Ist diese Bedingung erfüllt, so kann, wie G. DAHLQUIST?) gezeigt hat, die Interpolations- 
formel (1) mit der TAytor-Entwicklung von y(&;+1) an der Stelle x —= x; höchstens bis zum 
Gliede ya) - hr+2 einschließlich übereinstimmen. 


Für solch eine Übereinstimmung benötigt man in (1) offenbar m + 3 Koeffizienten. Ins- 
gesamt stehen 2m + 3 Koeffizienten in (1) zur Verfügung. Dain einer Interpolationsformel das 


; , , 2m+2 
Glied Ayryı Yarı Stets auftreten muß, kann man also | on Interpolationsformeln vom 
m 
Approximationsgrad h”+2 aufstellen. 


. Von diesen Interpolationsformeln interessieren natürlich nur solche, die die obige Stabilitäts- 
bedingung erfüllen. 


3. Fehlerfortpflanzung der Interpolationsformeln 


Für die praktische Anwendung sind numerisch stabile Interpolationsformeln mit möglichst 
günstiger Fehlerfortpflanzung von Interesse. Wir wollen daher versuchen, einen die Fehlerfort- 
pflanzung charakterisierenden Parameter zu finden, der es uns erlaubt, aus der Zahl von numerisch 
stabilen Interpolationsformeln solche möglichst günstiger Fehlerfortpflanzung auszuwählen. 


Bezeichnen wir den wahren Wert der Lösung unserer Differentialgleichung erster Ordnung 
an der Stelle x = %;—. mit y(x—,) und den Wert, den unsere Interpolationsformel (1) an dieser 
Stelle ergibt, mit ys—,, so erhalten wir aus (1) für den Fehler Y—u — Ya.) = &—„ bei Ver- 
nachlässigung der mit h multiplizierten Summe: i i 


m 


&+1 = L Ar u u—T AT N A Er (4), 
ins re rs 


!) H. RurisHAuser, Zeitschr. £. angew. Mathematik u. Physik III (1952), S. 65— 74. 
2) G. DAHLQUısT, Math. Scand. 4 (1956), 8. 33—53, 
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wobei: 
3 
T=C-ya)- hm+3 EEE EIER RO REHE GERNE 


das erste Glied des Abbruchfehlers unserer Interpolationsformel bedeuten soll. 
Macht man für den Fehler e, einen zu (5) analogen Ansatz: 


En lm Amis a (6) 
so folgt aus (4): 
m 
rn ZA —C RE TR. Sö5, Er een fo (7). 
u= 


Wir betrachten zunächst die homogene Differenzengleichung (7). Der Lösungsansatz 
&%—„ —*=# führt (nach Multiplikation mit A”—%) direkt auf die Gleichung (3), die aber für die 
uns allein interessierenden numerisch stabilen Interpolationsformeln (1) nur Nullstellen A mit 
[Al < 1 besitzt. Wegen |A| < 1 ist die allgemeine Lösung der homogenen Differenzengleichung (7) 
für die Fehlerfortpflanzung von nur geringem Interesse, und die noch hinzuzufügende partikuläre 
Lösung der inhomogenen Gleichung (7) ist für die Fehlerfortpflanzung ausschlaggebend. 


Der Ansatz: 
SE SKATER DEREN, GN DE NOBERR (5) 


für diese partikuläre Lösung führt unter Berücksichtigung der TAyLor-Entwicklung von (1) zu: 


ee (1 + &u- Ar.) RE NE NER re 


Die von den Abbruchfehlern der Interpolationsformeln (1) herrührenden resultierenden 
Fehler werden also nach einer größeren Schrittzahl näherungsweise im Verhältnis der Parameter E 
dieser Formeln stehen. 


Diese eigentlich nur für h = 0 zutreffende Aussage wird indessen im allgemeinen auch noch 
für kleine h + 0 angenähert richtig sein. 
In E haben wir damit einen praktisch leicht zu ermittelnden Parameter gefunden, der uns 


gestattet, beliebige numerisch stabile Interpolationsformeln z.B. mit den klassischen Interpola- 
tionsformeln von J. C. Anams in bezug auf Fehlerfortpflanzung zu vergleichen. 


4. Numerisch stabile Interpolationsformeln mit günstiger Fehlerfortpflanzung 
für m=1bis m=6 


2m+2 . r 
Die oben erwähnten | i ns) Interpolationsformeln vom Approximationsgrad h”*+? sind 
m 


systematisch in bezug auf Stabilität und Fehlerfortpflanzung untersucht worden. Nach Berech- 
nung der Koeffizienten A;,_„und A;_„ wurden die Interpolationsformeln nach (3) auf Stabilität 
geprüft und für die stabilen Interpolationsformeln der Fehlerfortpflanzungsparameter E be- 
rechnet. 

Die folgenden Tabellen 1 und 1’ geben die Koeffizienten A, „bzw. A;_ „für die numerisch 
stabilen Interpolationsformeln mit kleinstem E-Wert für m = 1 bis m = 6 wieder. 


Tabellel. Ar—u 


m Ar Ar—ı Ak—2 Ak—3 Ar—4 Ak—5 Ak—6 
4 l 
: 5 r3 
- 9 9 en 
1% 17 17 
16 il 
13300 8775 756 | 
4 0 v 21319 21319 21319 
243 1 79 
Zr = 0 0 a 
1000 . 8 125 
28107625 N 3568131 434875 1994 625 
0 x 82490048 10311256 1288907 82490048 


“. E i ” ur ee a rn BE PETE RT, Pe 
> ; - u Ku En es 


% wi 
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Tabelle I’. Ay 


18 0 
17 
1 10 2 2 
s 3 ry e 3 9 
6720 29700 13500 21300 R a 
i 31319 21319 21319 31319 
3 567 3 81 9 
u ä B, 0 = buch 
s 10 400 9 3 i 80 50 
379680 | 63007875 | 330750 | 34489815 2 35317485 t 5 
® 1288007 | 21945024 | T288007 | RR 41245094 


Das Verhältnis der E-Werte dieser Interpolationsformeln zu den E-Werten der Apamsschen 
Formeln ist in Tabelle 3 wiedergegeben. 


Bei der Auswahl dieser Interpolationsformeln wurden auch solche Formeln unberücksichtigt 
2 7 - i 
gelassen, die wohl für h = Osstabil sind, jedoch schon für kleine Werte - -h= O instabil werden. 
Ob dies eintritt oder nicht, hängt davon ab, ob die zu eingeschleppten Lösungen von (1) gehörigen 
Nullstellen des Polynoms 


M 
3 (v4 Ai) ar =0 ee 
u=—l1l y 
für kleines h= 0 aus dem Einheitskreis herauswandern oder nicht. 


Bei der Aufstellung der Tabelle 1 und 1’ sind nur solche Interpolationsformeln zugelassen 


worden, deren eingeschleppte Lösungen zu Nullstellen von (10) führen, die für — hi <0,1 


sämtlich absolut genommen noch kleiner als 0,95 sind. Der Zahlenwert 0,95 ist natürlich etwas 
Willkürliches. Wir wollten damit nur ausschließen, daß unsere Interpolationsformeln zu Unregel- 
mäßigkeiten in den Lösungswerten führen können. Denn solche Unregelmäßigkeiten können, 
wenn die zu eingeschleppten Lösungen gehörigen Nullstellen von (10) absolut genommen nur 
wenig unterhalb 1 bleiben, bekanntlich unter Umständen auftreten. 


In Tabelle 2 geben wir für die Interpolationsformeln der Tabelle 1 und 1’ den Absolutbetrag 
der absolut genommen größten zu einer eingeschleppten Lösung gehörigen Nullstelle A für 


[%) 
Zur h=-+0,1an 
öy 
Tabelle 2. |Almax 
a! m—=2 m 8 m—=4 m=d | 1 
öf r 
Mr h=—0,L1 0,2125 0,6003 0,7704 0,8838 0,9281 | 0,8725 
[7] 
a h=-+0,l1 0,1885 0,5474 0,7736 0,8985 0,9390 0,9276 


5. Beispiele 


„Um zu zeigen, daß auch noch für kleine h + 0 die Fehlerfortpflanzung durch den Ausdruck 
(9) für den Parameter E gut charakterisiert wird, haben wir als Beispiele die Differentialgleichung 


Y=-j-y (>) ee 


mit der Anfangsbedingung: x, = 0, y (Lösung: y = e*) 


u u TED ERBEN! 


und die Differentialgleichung 


a FE Eee Re Te 
B Re men BR R ü ber MD f LT > „N 2 4; 
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Y-/-—_ <0). ee 


mit der Anfangsbedingung: o=-l,y=1 (Eösungzs = l|x) 


nach Apams sowie nach den Interpolationsformeln der Tabelle 1 und 1’ numerisch integriert. 


In Tabelle 3 sind die theoretischen Werte E/E,nams mit den Ergebnissen der numerischen 
Integration nach 50, 150 und 250 Integrationsschritten verglichen. Als Schrittweite der Inte- 
gration wurde h = 0,01 gewählt; die mit einem Stern (*) markierten Werte wurden jedoch aus 
Integrationen mit der Schrittweite A = 0,02 ermittelt, da die Schrittweite h = 0,01 in diesen 
Fällen zu Abbruchfehlern führten, die bereits ebenso klein waren wie die Abrundungsfehler der 
Maschine?). 


Tabelle 3. Z/Eanams 


Aus numerischer Integration 


m Diff-Gl. He nach N Schritten 

N=50 N=150 N = 250 
| en | ns | omese 
|| 0 | ae an | sie 
>| | | | oa | oa 
|| | en | 0220 | 0220 
s| a | an | om | Doros | .oloros 
| | me | oa | one | one 


Die Abweichungen der durch numerische Integration gefundenen Fehlerverhältnisse von den 
theoretischen für A = O gültigen Werten betragen in unseren Beispielen für N = 250 und h = 0,01 
weniger als 1%. Die theoretischen Werte (9) bieten also einen auch für praktische Zwecke brauch- 
baren Maßstab für die Auswahl von Interpolationsformeln (1) mit günstiger Fehlerfortpflanzung. 

Mit zunehmendem m gelangen wir in Tabelle 3 zu immer günstigeren Formeln bezüglich 
Fehlerfortpflanzung, bis sich schließlich für m = 5 für unsere Interpolationsformel (Tabelle 1 
und 1’) der Fehler auf nur noch 8%, des Fehlers der entsprechenden Anams-Formel reduziert. 


Daß für m = 6 der günstigste Wert E/Eynanms wieder anwächst, hängt wohl damit zusamm- 
men, daß mit zunehmendem m immer mehr eingeschleppte Lösungen der Differenzengleichung (1) 


. auftreten und sich dadurch die Tendenz zur Instabilität vergrößert. 


Teil 2: Formeln für Differentialgleichungen zweiter Ordnung 


6. Numerische Stabilität der Interpolationsformeln 


Wir übertragen die Stabilitätsuntersuchungen von H. RuTISHAUSER!) auf unsere Formeln (2). 
Wegen des Auftretens der Glieder mit y;—„in der ersten Formel (2) können wir die Ergebnisse 
von RUTISHAUSER nicht direkt übernehmen, sondern müssen sie auf unseren Fall ausdehnen. 

In RUTISHAUSERS Schreibweise erhalten wir aus unserer Differentialgleichung zweiter Ord- 
nung y” = f(x, y, y') und aus unseren Interpolationsformeln (2) die Variationsgleichungen: 


of a E 
N ———. Nr _— 1 — 0 R 
ER Net ER ML NE 
1 m Up A nn u 
74 x > Ay NKk—u a en NRk—u + h- 2 Ara Nk—u — 0 ’ Umso  ä (13) 
u=—l De BE 


m m r ’ 6 2 177 r 
> Bi, Nk—u zz h- ER De, Nk—u + h2. > DB; Ne—u = 0 
u=—1 


u=—l u=—l1 


wobei Aryı =— 1; Ber ı = — 1und Bi4ı = 0 zu setzen ist. 


3) Die Rechnungen wurden auf einer IBM-704 Maschine mit einer 20 Dezimalen liefernden interpreta- 
tiven Subroutine des National Bureau of Standards, Washington, D. C., durchgeführt. 
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Der Ansatz: 


NE—a = Po Are, Ma = Pı- Ar, 9 (14) 
führt zur folgenden Determinantenbedingung für A: 
| 4 U = 
&y ay. 
RE 1 Sa, ZA ,Mul=0...(). 
h BE i u Re u=—l cp Fr l 


m m - 
SB.4m 0 h- SZ BoAo My Bl,“ 
um—L we uu—1 
Multipliziert man die beiden letzten Zeilen der Determinante mit 4”—* und entwickelt die 
Determinante, so erhält man, wenn man nach Potenzen von h ordnet: 


m m m m 
> ER nd > Be, AP — P3 N PR aa oe > Di 


u=—1 u=—1 n=—1 ne=—1 
öf m m m ur 2 

+I.n.]) 2 ame ID BA m— 3 Ay, Mma. DZ Biel Hs —=0lIß): 
oy' u=—1 u=—1 u=—l u=—l1l 


Aus (16) folgt dann für genügend kleines A als hinreichende Stabilitätsbedingung: 
Die Interpolationsformeln (2) sind numerisch stabil, wenn das Polynom 


SZ Ay .m=r. 5 B_,imme— 3 Adam. 5 B_,Ie#=0...(1 
1 


u—rZl u! ie n=—1 
außer der trivialen doppelten Nullstelle A = -+ 1 nur Nullstellen A mit |A| < 1 besitzt. 


Unser Kriterium (17) ist für die numerische Auswertung komplizierter als das entsprechende 
Kriterium von RUTISHAUSER. Dies hat seinen Grund darin, daß letzterer die Interpolationsfor- 
meln für y und für y’ einzeln auf Stabilität hin untersuchen konnte. Durch Einführung der Glieder | 
mit ye—„in die erste Formel (2) haben wir indessen unsere Formeln so gekoppelt, daß dies nicht 
mehr möglich ist. Der Vorteil dieser Kopplung wird jedoch aus Nr. 8 ersichtlich. 

j 
1 


7. Fehlerfortpflanzung der Interpolationsformeln 


Wir beschäftigen uns wieder mit dem Problem, wie sich die Abbruchfehler unserer Inter- 
polationsformeln im Laufe der Rechnung fortpflanzen, und wir suchen wieder einen Maßstab zu 
finden, der es ermöglicht, verschiedene numerisch stabile Interpolationsformeln hinsichtlich 
ihrer Fehlerfortpflanzung zu vergleichen und Formeln mit möglichst günstigerFehlerfortpflanzung 
ausfindig zu machen. 


Bezeichnen wir die wahren Werte der Lösung unserer Differentialgleichung mit y(;_.) 
und y’(@&—,„) und die Näherungswerte unserer Interpolationsformeln (2) mit yr_,„ und Yrp50 
ergeben sich für die Fehler: 


Y—u Ya) — Nk—u > Ye—u Ze Y(Rr— u) Fr Ne— u A ne (18) 


aus (2) bei Vernachlässigung der Glieder mit der höchsten A-Potenz die folgenden Differenzen- 
gleichungen: 


u=—l 


. 1 ” m R a 
N+1 NR R N Ar N—u + 2 Ara Nk—u —” Tr ’ | 
p= 


m m ee 7 a (19) 
Ne+r = 2, Bi. Mu th: ZB Na Tr | 
um u=0 


“mr m . . r 
wobei T; und T;, die jeweiligen Abbruchfehler unserer Interpolationsformeln für y’ und y an der 
Stelle x = x; bedeuten sollen. 


Durch Umordnung von (19) erhalten wir, wenn wir noch 


Ben hm dm ee ee 
m m 
Nk+1 2 Bau NE—u & Bi—u —u TEE Tr > 
m s 2 1). 


dar A + Artı Be) Mu Ar + Aryı Br.) «u = — Arıı r—hTi 
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Unsere Interpolationsformeln (2) mögen y bis zum Gliede h? der Tavıon-Reihe und y’ bis 
zum Gliede AP—1 der Tayror-Reihe genau erfassen. Dann gilt für unsere Abbruchfehler T, und 
T; in erster Näherung: 


T=C- ya) bett, IE AN ne (22) 


wobei C und €’ Konstanten sind, die jedoch von Interpolationsformel zu Interpolationsformel 
verschieden sein werden. 


Macht man einen (22) entsprechenden Ansatz: 


ut, ehe. nn... (23) 
so folgt aus (21): £ 
m m 
ED 2 De kR—=—C, 
n=0 u=0 
m in (24). 
Pr+1 Ara + Aryı Br_,) 5. — (Ar + Arzı Br.) * Bra = — Aryı CC’ 
= Pe 


Wir betrachten zunächst wieder die homogenen Differenzengleichungen (24). Machen wir 
für diese homogenen Differenzengleichungen den Lösungsansatz: 


Kk—u =0&° ART P Pen = ß a N N Eh: (25) 
so finden wir für } wieder unsere Gleichung (17) und können wie bei den Interpolationsformeln für 
Differentialgleichungen erster Ordnung wieder schließen, daß für die Fehlerfortpflanzung nume- 


risch stabiler Interpolationsformeln die partikuläre Lösung der inhomogenen Gleichungen (24) 
ausschlaggebend sein wird. 


Für diese partikuläre Lösung machen wir den Ansatz: - 
5 =D-k+E-R, PB ST a a 26): 
Setzen wir (26) in (24) ein, so erhalten wir durch einen Koeffizientenvergleich der Glieder 
mit k2, der Glieder mit k und der freien Glieder Bedingungsgleichungen für die Koeffizienten D, 
E,F. 
Unter Benutzung der TAYLorR-Entwicklung von (2) sieht man leicht, daß die Bedingungs- 


gleichungen für die Glieder mit k? identisch erfüllt sind. Eine der beiden Bedingungsgleichungen 
für die Glieder mit k ist ebenfalls identisch erfüllt, die andere führt zu: h 


EA re ee (27). 


Aus den beiden Bedingungsgleichungen für die freien Glieder kann man dann die verbleiben- 
den Koeffizienten berechnen. Die elementare Rechnung ergibt für den für das Wachstum von a; 
ausschlaggebenden Koeffizienten E den Ausdruck: 


(An — Eu: A) +C(1+ Zu: Bi-) 
u= u= 


Ro). 


I se m. , m R = m Mm PR 
(Ars — Zur): 3 Bit ( I Bu): > Ai 
wa 


v=1 B=— al 


8. Numerisch stabile Interpolationsformeln mit günstiger Fehlerfortpflanzung 
fürm=1bis m=3 


Wir haben für m=1 bis m = 3 wieder eine vollständige Untersuchung aller möglichen 
Kombinationen von Interpolationsformeln (2) für Differentialgleichungen zweiter Ordnung 
bezüglich Stabilität und Fehlerfortpflanzung durchgeführt mit dem Ergebnis, daß wir noch für 
p=m-+4 numerisch stabile Formeln (2) erhielten. Demgegenüber besitzt die Anamssche 
Interpolationsformel für Differentialgleichungen zweiter Ordnung nur einen Genauigkeitsgrad 
p=m-3. Sich auf gleichviel rückwärtige Stellen stützend, gewinnt man also gegenüber 
Anpams eine h-Potenz an Genauigkeit. 

Unter der Annahme, daß die interpolatorischen Glieder Ay+ı yg+ı und Br+ı Yr+1 in (2) 


smbMN A ee 3m-+ 
stets auftreten, bestehen offenbar | “ Möglichkeiten für die erste Formel (2) und | “ 
m + 


en 


mögliche Kombi- 
m-4 m-4 


Möglichkeiten für die zweite Formel (2), d. h. insgesamt | 
nationen (2). 

Für m = 1 sind dies 126, für m = 2 bereits 17640 und für m=3 schließlich 1359072 Kom- 
binationen. 


rim - 
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Die folgenden Tabellen 4, 4 und 4” geben die Koeffizienten für die numerisch stabilen 
Interpolationsformeln (2) mit kleinstem E-Wert (28) für m = 1 bis m = 2 wieder. 


Tabelle 4. Ak—ı Br—u 


m Ar+ı | Ar | Ar—ı | Ak—2 Ak—3 | Br | Bk—ı | Bı—2 | Bk—3 
23 4 7 r - 
: IB = IB 
11 9 8 - 157 m 
2 6 0 Be} 3 66 66 0 
1352 2376 | 1024 r | 16000 ® 704 227 
2 1007 . 1007 1007 9891 1099 9891 


Tabelle 4. Ak—n, Bt—u 


m Ar Ar Ak—2 Ak—3 | Bk Bı—ı Bı—2 Bk—3 
I  —— 
4 1 | 8 14 
1 == er 15 15 
| 98 2 1 
2 0 0 ) | 32 an 
| 
3 ) 0 az Sl | 2176 m | go 0 
1007 1007 3297 1099 


Tabelle 4". Ak, Bk—u 


N RE ER een VE | Beh | B | Bi-ı | Bio | Bi-ı 
| | 
1 | 4 1,8 Io 
en £ + a 
| I j 
ER! 2 a Te 239 „| 3 |;., 433 
5 Be 5 5 | 1188 | 1086 | 9504 
8 222 004 204372 r 126 68 | 1812 | 956 0 
1007 1007 | 1007. | 1098 | 1099 , 109 


Über die Wanderung der A-Werte der eingeschleppten Lösungen von (2) für A =+ 0 entschei- 
den jetzt die Nullstellen von (16). 

Bei der Auswahl der Interpolationsformeln (Tabelle 4, 4°, 4”) wurden die gleichen beschrän- 
kenden Bedingungen gestellt wie im Falle der Differentialgleichung erster Ordnung, nur daß 


7 - han die Stelle von er -htritt. 
oy oy 


In Tabelle 5 geben wir die der Tabelle 2 entsprechenden Werte |Almax für die obigen Inter- 
polationsformeln (Tabelle 4, 4, 4°’) wieder. 


jetzt natürlich 


Tabelle 5. |Almax 


m=|1l m=%2 m=3 
a 
of | 
a h=—0l 0,8751 | 0,7607 0,9311 
of | 
y' h=+01 0,8017 | 0,7615 0,8985 


9. Beispiele 


Wir bringen wieder drei Beispiele, um zu zeigen, daß der Ausdruck (28) für E auch fü 
h=# 0 die Fehlerfortpflanzung gut charakterisiert: n ” a 


[77 ’ ö 
a N (0) a 


E. Fratsere, Numerisch stabile Interpolationsformeln mit günstiger Fehlerfortpflanzung 109 


mit der Anfangsbedingung: u =0, y=1,y=—1 (Lösung: y = e®), 


De 1 of = 
MB, Re en 
mit der Anfangsbedingung: u = 1, =1, =, (Lösung: y = Vz), 
7 y? öf 
—f=-2.3- ee: 
| He : E (31) 
mit der Anfangsbedingung 9 = 1, =1L,,=—1 (Lösung: y— a 
2 


In Tabelle 6 sind die theoretischen Werte E/Eynıms und die Ergebnisse der numerischen 
Integration (Schrittweite h = 0,01) für diese drei Beispiele sowohl für die Fehler in y als auch 
für die Fehler in y’ zusammengestellt. 


Tabelle 6. Z/Eananms 


Diff. | Theoret. Wert Aus numerischer Integration nach N Schritten 


eier nach (28) für y für y’ 
N=50|N=10|N=-20| N=590|N=10| N=-250 


29) 0,2409 0,2386 0,2384 | 0,2390 0,2334 0,2384 
30) 0,2368 | 0,2367 0,2345 0,2342 0,2342 0,2338 0,2338 
3l 


( 
ı I( 
(31) 0,2338 0,2318 0,2315 0,2296 ' 0,2047 0,2328 


(29) 0,0074 0,0050 0,0048 0,0094 0,0050 0,0048 
2 | (30) —0,0011 0,0013 |—0,0008 |—0.0012 0,0004 —0,0015 |—0,0016 
(31) —0,0144 |—0,0157 —0,0158 | —0,0144 —0,0021 |—0,0154 


(29) —0,0061 |—0,0062 —0,0063 | —0,0033 —0,0062 |—0,0063 
3 (30) 0,0055 0,0038 0,0047 0,0049 0,0054 0,0051 0,0053 
(31) —0,0038 —0,0027 |—0,0025 | —0,0002 —0,0127 —0,0033 


Die in Tabelle 6 zum Vergleich herangezogene Apamssche Formel stützt sich auf eine rück- 
wärtige Stelle mehr als unsere Interpolationsformeln. Denn erst dann erreicht man mit den 
Formeln von Apams den gleichen Genauigkeitsgrad wie mit unseren Formeln. 

Bemerkenswert ist der geringe Fehler unserer Formeln fürm =2 und m=3; er beträgt 
für h= 0 nur etwa 19/,, bzw. 5°/,. des Fehlers von ADAms, und seine Größe wird bei der nume- 
rischen Integration von (30) (& h= 0) gut bestätigt. Für die Beispiele (29) und (1) | en h 0) 
treten dagegen wie bei den Formeln für die Differentialgleichung erster Ordnung wieder Ab- 
weichungen von den theoretischen Fehlerwerten von der Größenordnung von 1% auf. 

Wir haben für m=1 und m=2 auch Interpolationsformeln (2) untersucht, bei denen 
in der ersten Formel die erste Summe von u = O bis u = m und in der zweiten Formel die zweite 
Summe von u = — 1 bis u = m erstreckt wird. Die sich dann ergebenden Formeln sind jedoch 
nicht so günstig wie die oben mitgeteilten. 


10. Der Spezialiall y” = f(x, y) 
In einer neueren Arbeit hat G. DAutquıst®) u. a. bewiesen, daß für die Differentialgleichung 
y’ = f(a, y) für ungerades m numerisch stabile Interpolationsformeln von der einfachen Form: 


m Mm 
Yr+1 —— Die YUr—u + h2 ° >) Bi YUr—u ar wohn EM ee ke Me (32) 
u=0 


u=—l 
existieren, die wieder den Genauigkeitsgrad p = m + 4 besitzen. 
In Tabelle 7 und 7” haben wir fürm = 1,m = 3und m = 5 die Koeffizienten der numerisch 
stabilen Interpolationsformeln (32) mit günstigster Fehlerfortpflanzung (x; = E- K2) angegeben. 


Tabelle 7. Bi — u 


m | Br | Bk—1 | Br—2 | Bk—3 | Bk—4 | Bk—5 
1 2 —l 

3 1 0 1 = 

5 1 0 0 0 1 —l 


#) Thesis: Stability and error bounds in the numerical integration of ordinary differential equations, 
Stockholm 1958. 
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Tabelle 7”. Bk—u 


Ba: | 
. 12 s 12 | 

17 29 111 9 | es 
3 240 30 10 FT 71) 

787 2027 3263 3751 | 3263 2027 187 * 
5 | 5006 | 6 4033 304 , 08 | RG | 706 


"Für m =1 ergibt sich die wohlbekannte Formel der zentralen Differenzen, die hier also 
im Gegensatz zu dem Fall der allgemeinen Differentialgleichung y” = f(&, y, y’) stets numerisch 
stabilist. Fürm = 1 ist dies die einzige Formel, die das Verlangte leistet. Für m = 3 und m = 5 
existieren dagegen mehrere solche Formeln, aus denen wir die angegebenen als diejenigen mit 
kleinstem Fehlerparameter E ausgewählt haben. Wir haben hier wieder an Beispielen bestätigt, 
daß das Verhältnis der E-Werte ein recht genauer Maßstab für die relative Fehlerfortpflanzung 
solcher Formeln ist. 
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Über periodische Lösungen nichtlinearer Differentialgleichungen 
und über eine Beziehung zwischen Amplitude und Frequenz 


Von P. SAGIROW 


Die Arbeit ist eine Fortsetzung von [3] und behandelt die Fehlerabschätzung bei der Berechnung der perio- 
dischen Lösungen der Differentialgleichung L(D)xz = f(x, &) mit Hilfe der harmonischen Balance. Für die 
Frequenzen beliebiger periodischer Lösungen von L(D)x = f(x) wird eine obere Schranke angegeben. 


The paper is a continuation of [3] and treats the error estimation for the calculation of the periodic solu- 
tions of the differential equation L(D)x = f(x, x) with the aid of the harmonic balance. An upper bound is 
given for the frequencies of arbitrary periodic solutions of L(D)x = fix). 


Pa6oTa ABIAETCH IPONOJBKEHHEM [3] H ITOCBAIMEHA OMCHKE OIIMÖRH IPu pacyere TIepHoNHu- 
yYeckux pemeHnä AudbhepeHnnmmaspHoro ypasHueHnun L(D)e = f{x,&) C HOMOMBIO TAPMOHHYe- 
CKOTO ÖaslaHca. YRkasaHa BepxHAA TPAHB IS YACTOT MOÖBIX IIePHONHYecKHX pemmeHnü ypaBHe- 
HuA L(D)e = f(&). 


In [3] wurde der Fehler abgeschätzt, mit dem die Näherungswerte der Frequenz und Ampli- 
tude sinusähnlicher periodischer Lösungen einer nichtlinearen Differentialgleichung behaftet sind, 
wenn man diese mit Hilfe der harmonischen Balance löst. Untersucht wurde der einfachste 
Fall L(D) x + f(&) = 0 mit einer eindeutigen und ungeraden Funktion f(&). 

In der vorliegenden Arbeit wird die Ausdehnung des in [3] vorgeschlagenen Verfahrens auf 
allgemeinere Fälle betrachtet. Eine solche Verallgemeinerung ist in zwei Richtungen möglich: 
a) Man kann kompliziertere nichtlineare Differentialgleichungen ins Auge fassen und 

b) man kann über die erste Näherung (sinusähnliche Lösungen) hinausgehen und höhere 
"Näherungen betrachten, bei denen neben der ersten auch die höheren Harmonischen berücksichtigt 
werden müssen. 

Die Frage a) wird am Beispiel der Gleichung L(D) x = f(x, &) im ersten Teil des Beitrages 
behandelt. Da die Untersuchung im wesentlichen nach dem in [3] vorgeschlagenem Schema ver- 
läuft, begnügen wir uns mit einer knappen Darstellung. Die Ergebnisse können ohne Schwierig- 
keiten auf den allgemeinen Fall KKD)x = R(D) f(x, &, &, ...) übertragen werden. 

Det zweite Teil ist der Frage b) gewidmet. Es wird u. a. gezeigt, daß die Frequenz der perio- 
dischen Lösung der nichtlinearen Differentialgleichung L(D) x = f(x) einer algebraischen Glei- 
chung genügt, deren Koeffizienten in einfacher Weise von den Amplituden der einzelnen Harmo- 
nischen abhängen. Eine Untersuchung dieser Gleichung zeigt, daß die Frequenzen aller (nicht: 
nur der sinusähnlichen !) periodischen Lösungen von L(D) x = f(x) stets kleiner sind als die größte 
positive Wurzel der Gleichung Y,(o) = Im L(io) = 0. 
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I. 
Betrachtet werde die Gleichung 
LUD)z=f@), WAD)=r,D®r +7, Dm-1il..4n,_1ı Dtm ee lb) 
wobei die partiellen Ableitungen von f(z, &) beschränkt sein sollen: 
fx, &) If@; &) 
— trEeM Ze Me (2). 


Sucht man eine periodische Näherungslösung x = = + asin » tmit Hilfe der harmonischen 


Balance, so erhält man mit unseren Bezeichnungen zur Bestimmung von o, a,und a das Näherungs- 
system [1] 


Img = Ay; aX=B,, ae TA Ener (3), 
wo 
27 2n 27 
— ji _ 1 — T 
A=,, e-- -) du, A, =— | /(- -)cosudu, 3-2 | e-Jsinwau 
7 7 7 
ö ö ö (4) 


fe = + asinu, @ acos u) 


ist, während X und Y Real- und Imaginärteil von L(io) sind: L(io) = X(®) +i Y(w). Die 
Fehler, mit denen die aus (3) erhaltenen Näherungswerte ®,, Qy,, a, behaftet sind, bleiben un- 
bekannt. 

Um diese Fehler zu bestimmen, verallgemeinern wir das in [3] vorgeschlagene Verfahren. 
Wir machen den Ansatz 


x) = + asin 917-3 (a, cosnot Fb,sinnel). ......Rh. (5) 


n=2 
(beim Ansatz in [3] war a, = 0 und es wurde nur über ungerade Indizes summiert) und ent- 
wickeln f(&,&) in eine Fourikrreihe 


ja.) = + BEA CO No tus, Enz ee En (6). 


n=1 
Die rechte Seite von (6) wird an Stelle von f(x, X) in (1) eingesetzt und die entstehende inhomogene 
lineare Differentialgleichung gelöst. Ein Koeffizientenvergleich zwischen der-erhaltenen Lösung 
und dem Ansatz (5) liefert zur Bestimmung von ©, 4y, d, An, D„ das unendliche Gleichungssystem 
7. = Ag ax, =D, aYyy=4, 
Ye = Dr; VE RE TTR DE, 
m &n + Dn Yn — A, 
wor. — (nl)... =. 2.0140); ist, 
Die Fourierkoeffizienten von (6) in den ersten drei Gleichungen (7) formen wir um, wobei 
wir der Kürze halber die Bezeichnungen benutzen 


[6.0] [0,0] 
u=ot, Z=N („cosnot+b,„sinneol), 2’=N(—- na, osinnoet—+nb„ wcosnol) 
2 n=2 2 n=2 
Dann ist z.B. 
27 


1 


2n 
„|t@+emu+z, awcosu+ Zr)du= 5, | I + asinu, a 0 cos u) du 
2 2 2 an, 2 
Ö 


A= 5; 
1 27 
I N [28 + 27 rn] du. 


Das erste Integral ist A, aus (4), das zweite Integral bezeichnen wir mit &, und schätzen es ab. 
Auf diese Weise erhalten wir 


A=A+%: AHA, +0; B=B+B; dr . (8) 
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wobei Ay Ay B, durch (4) gegeben sind, während 


ä 4 ri 
|xo]| £ M,ı max 2 + M; max |}|; ll, Ifı] = — Mı max 2 + M, max |Y,) 
2 


ist. Die Beträge | 2| und | 2 werden ähnlich wie in [3] abgeschätzt. Es ist 
2 2 
Br SA Be 
It < — 3 —— - V fl(elu), Hu 
zIszja+=2, 7) Sm Zur, Ka), &W) 


wo V die totale Variation auf dem Intervall [0, 2 =] der u-Achse bedeutet und L, = |Un i o)| ist. 
Berücksichtigt man die Ungleichung V f(x(u), 2(u)) s M,V x(u) + M, V (u) und beschränkt 
man die Betrachtung auf diejenigen o-Werte, die der Bedingung 


52 (m,2 I 2m, 0 3 As IE (9) 


n=24n n=2 L 


genügen (vgl. die entsprechende Bedingung (13) in [3]), so erhält man auf demselben Wege wie in 
[3] endgültig die Abschätzungen 


1 u, +0 M) 


= lu ne me Se Sg Fe (10) 
2 = Sr 2 A or ") S,nL, 
Be 12, In £ n=2 L, 
und 
4 
\n|<x*a; ll, |ßıl = 0, 2 je ee ee (11), 
wo zur Abkürzung | 
8 y2 
—M+Ma) .SM+Mo (12) 
8 y2 n=2 n Lu e 


St © S 
ee a, Bu! 
7 (IM Z, In an E > L, 


gesetzt ist. 


Damit sind alle Abschätzungen gewonnen, die zur Anwendung des in [3] geschilderten 
Verfahrens auf die Gleichung (1) benötigt werden. Die Bestimmung der Intervalle, in denen die 
strengen Werte von ®, ay, a liegen müssen, geschieht unter Berücksichtigung der Abschätzungen 
(11), (12) nnd der Darstellung (8) aus den drei ersten Gleichungen des Systens (7) 


Mo=At%: aN=B +ßı and iu Sera 


Da (13) aus drei gekoppelten transzendenten Gleichungen besteht, ist die Lösung im allgemeinen 
mit sehr großem Rechenaufwand verbunden. Das trifft allerdings auch auf das in der Theorie 
der automatischen Regelungen benutzte Näherungssystem (3) zu, wo. = & = ßı = 0 gesetzt ist. 

Wir weisen auf drei wichtige Sonderfälle hin, in denen sich die Rechnung wesentlich verein- 
facht. 


2an 
E x — 1 s 
Falll /(&,2)=/f,(). Dann ist A, = a Ih Fr + asin u) cos udu=0 und deshalb 
ö 


_ 


4 ne 4 | . 
| Yı(o)| < 7% Der Spezialfall, wo f(x) außerdem ungeradeist, wurde in [3] ausführlich behandelt. 


Fall II: f(x, &) = f,(&). Dann ist B, = O0 und folglich ist |X,(o)| < ar 3 
1 


Fall III: f(x) = fs(Ax + u&); A, u = const. In diesem Fall ist 

2n 

Fan 1 i 

A, = = | Ta A % + asin u) + ua wcos ı cos u du 
(1) 
27 

_ 1 Ay n : 

B=-|f \ & + asin u) + ua wcos ı sin u du 


TT 
0 
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und deshalb Aa A, —uaoB,—=0. Dann ist aber 
4 
A Yılo) — u 0 Xo)| <x(A| + | 0). 


In allen drei Fällen kann ® unabhängig von a, und a bestimmt werden (vgl. [3). An- 
schließend können a, und a aus zwei Gleichungen graphisch ermittelt werden. | 


I. 


Die harmonische Balance dient in erster Linie zur Bestimmung sinusähnlicher periodischer 
Lösungen. Es ist aber schon mehrfach versucht worden, mit diesem Verfahren auch allgemeinere 
periodische Lösungen zu ermitteln [2]. Es entsteht die Frage, wie man in einem solchen Fall eine 
Fehlerabschätznng erzielen kann. 

Der Ansatz (5) bleibt derselbe. Man entwickelt f(x, &) in die Fourikrreihe (6), löst dann (1) 
und erhält durch einen Koeffizientenvergleich das System (7). Will man außer der ersten Har- 
monischen k weitere Harmonische berücksichtigen, so muß man in (7) die ersten 2 k+3 Gleichun- 
gen betrachten. Die Aufspaltung der Fourıkrkoeffizienten in einen strengen und einen abge- 
schätzten Teil geschieht wie oben. Die praktische Bedeutung des erhaltenen Systems ist jedoch 
gering, da die Lösung auf unüberwindliche rechnerische Schwierigkeiten stößt. Wir verzichten 
deshalb auf nähere Angaben und benutzen das allgemeine System (7) lediglich zur Ableitung einer 
Beziehung, die uns neben einem wertvollen Einblick in die Struktur periodischer Lösungen eine 
praktisch wichtige obere Schranke für die Frequenz dieser Lösungen geben wird. 

Gegeben sei die Gleichung 


DIA ASCHHT ER ERTL ART (14) 


mit der periodischen Löusng (5). Die Größen ®, a, 4, An, d„ genügen dem System (7), wo 


2n 27 
An= | 1(@ + asinu + Z)cosnudu, Bu [1(@+asinu+ 5) sinnudu 
IT 2 2 
ö ö 
: AR 
ist. Es sei a f@). 


Wir bilden den Ausdruck 


zn 
oo 1 , i 
aA, +3(- na, B„+nb,An) = [1(@+ asinu + F)(acos u + z/)du 
n=2 2 P} 
(0) 
an % Zasinu 
e* [ . du = .n 0‘ 
TEN du TC 0 
(0) 
Es ist also a A, + 5 (— na, B„ + nb„A,) = 0. Setzen wir hier laut (7) A„ = nXn + Dn Ya» 
Bn = — % Yn + bu Xu ein, so erhalten wira Y, + Y n(an + b5) Yna=0 oder mit der Be- 
zeichnung c = @ + b„ die Gleichung 
Beer en Bello) 


n=2 


Diese können wir noch weiter umformen, indem wir Y, durch die Koeffizienten von (14) 
ausdrücken: Y, = m—_ın © — Im_3 (n ©)? + +». Dann bekommen wir nach einer Umordnung 


der Glieder die Gleichung 
(@ +22 +3? 0) Im—ı1® (a? 24 % 34 +) Tm—3 a Bi (16). 
a BE EEE 


ab- 


Man beachte, daß alle Koeffizienten von (16) konvergente Reihen sind, da die c„ mit 


nehmen (siehe oben Teil ]). a 1 
en Gleichungen (15) und (16) lassen eine wichtige Schlußfolgerung zu. Im Falle 


m = 3 erhält man aus (16) 


ul, 


p) 2 p2 2 22 ... h 
2+25+3°G5+ ee ER NE EL): 
a? 24 es 34 & en To To 


8 
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d. h. die Frequenzen aller überhaupt möglichen (nicht nur sinusähnlichen !) periodischen Lösungen 


r : 
von (14) mit m=3 sind kleiner als der Näherungswert @, = I den das Verfahren der harmoni- 


0 
hen Balance liefert. Ebenso verhält es sich bei Gleichungen vierter Ordnung. RS 
3 Dieses Ergebnis gestattet es u. a., die Fehlerabschätzung für das in [3] betrachtete Beispiel 
wesentlich zu verbessern. Dort wurde festgestellt, daß die Frequenz einer sinusähnlichen periodi- 


schen Lösung von &+2&+3& +42 +) = 0, |f@)| = 3 im Bereich y3 — 0,03 = o< 
Y3 + 0,03 liegen muß. Jetzt können 
wir behaupten: Die Frequenz einer 
sinusähnlichen periodischen Lösung 
der im Beispiel betrachteten Diffe- 
rentialgleichung liegt im Intervall 
Y3 — 0,03 < » < y3. Für die Am- 
plitude folgt dann im allgemeinen 
a > 4 

Aus (17) ist auch klar ersicht- 
lich, in welcher Weise das Auftreten 
von höheren Harmonischen die Fre- 
quenz beeinflußt. 

Um bei einer Differentialgleichung beliebiger Ordnung zu einer ähnlichen Aussage über die 
Frequenz zu kommen, betrachten wir die Gleichung (15). Denkt man sich L(i o) in der X+i Y- 
Ebene aufgetragen, so ergibt sich (falls L(D) ein Hurwırz-Polynom ist) eine spiralförmige Kurve 
[2], die den Ursprung umläuft (das Bild zeigt schematisch die Fälle mit m=8 undm=9). 

Die Schnittpunkte mit der X-Achse entsprechen den Näherungswerten @77, @n9 nz.» 
der harmonischen Balance. Ist & > @®,max, So sind alle Y,, Yz Y,... vom gleichen Vorzeichen, 
und (15) kann nicht erfüllt sein. Folglich sind die Frequenzen sämtlicher periodischen Lösungen von 
(14), wo L(D) ein Hurwitzpolynom ist, stets kleiner als die größte positive Wurzel &,max der Glei- 
chung Y(o®) = Im L(io) =. 


m=9 
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Thermokonvektive Zellularströmungen 
bei inkonstanter Erwärmung der Grundfläche 


Von JÜRGEN ZIEREP*) 
Herrn Prof. Dr. W. SCHMEIDLER zum 70. Geburtstag gewidmet 


Die Zellularkonvektionsströmungen über gleichmäßig erwärmten Flächen sind weitgehend untersucht 
worden. In der vorliegenden Arbeit wird das Problem für ungleichmäßige Erwärmung der Bodenfläche 
formuliert und für einen speziellen charakteristischen Fall sowohl im zweidimensionalen als auch im rotations- 
symmetrischen Fall explizit gelöst. Die berechneten Stromfelder sind schon wiederholt im Experiment 
beobachtet worden. 


The cellular convection currents over uniformly heated surfaces have been investigated to a large extent. 
In the present paper the corresponding problem for non-uniform heating is formulated and is solved for @ 
certain typical case both in two-dimensions and in the case of rotational symmetry. The computed flow fields 
have been repeatedly observed in experimental tests. 


Hesnuonstpupre KOHBERIMOHHBIE TIIOTORH Han PABHOMEPHO PA30TpeBAaeMbIMH HIOBEPXHOCTAMHU 
Y’Ke TTONPOOHO usyyeHst. B HacTosımeü pa6oTe aTa zanaya POPMyJIUpyercaH mA NOBEPXHOC- 
TeH, IONBePFAeMbIX HePABHOMEPHOMy Pa30TpeBy, U PemIaeTca IA HeKOTOPOTO XAPARTepHCTH- 
TECKOFO CJIyYAsI KAK B IBYXMePHOM TAK U BO BPAImaTeıbHo CHMMETPHYHOM cııyyae. Pemenung 
UOIYYAIOTCA B ABHOM Bulle. PacyerHsie HONA MOTORA HEeOAHOKPATHO HA6JONANHCL B OTBITE. 


*) Aus dem Institut für theoretische Gasdynamik der Deutschen Versuchsanstalt für Luftfahrt. 
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1. Einleitung und Übersicht 
; Die von Bänarp [1] entdeckten Zellularströmungen sind in den letzten Jahrzehnten 
wiederholt von der Seite der experimentellen und der theoretischen Forschung behandelt worden. 
Diese BEnarpschen Zellen bilden sich in einem zähen, wärmeleitenden Medium, wenn der Tem- 
peraturabfall von der Grund- zur Deckfläche einen kritischen Wert überschreitet. Von der Seite 
der Theorie wurde bisher ausschließlich der Fall einer gleichmäßigen Erwärmung der Grund- 
fläche (und konstanter Temperatur der Deckfläche) behandelt. Selbst bei diesen verhältnismäßig 
einfachen Voraussetzungen für die Theorie waren in der Darstellung von RAyLeic# [2] und 


 JEFFREYS [3] einige Fragen von grundsätzlicher Bedeutung offen geblieben, von denen einige 


z.B. durch eigene Arbeiten [4] in den letzten Jahren beantwortet werden konnten. (Bezüglich ein- 


‚schlägiger experimenteller Arbeiten hierzu siehe etwa [5].) Hier solles sich nun um eine Erweite- 


rung der Theorie insofern handeln, als wir das Phänomen der Zellularströmungen bei ungleich- 
mäßiger Erwärmung der Grundfläche (aber konstanter Temperatur der Deckfläche) studieren 
wollen. Diese Behandlung scheint mir für viele Anwendungen von ganz grundsätzlicher Bedeutung 
zu sein. Denkt man etwa an ähnliche Strömungsvorgänge in der Meteorologie, so stellt der Fall ° 
der inkonstanten Erwärmung der Grundfläche die Regel, der Fall der gleichmäßigen Erwärmung 
der Bodenfläche einen stark idealisierten Sonderfall dar. 


“ Was die Voraussetzungen unserer Theorie betrifft, so beschäftigen wir uns hier ausschließlich 
mit stationären Strömungen, obwohl wir die Grundgleichungen auch für den instationären Fall 
angeben werden. Wir schreiben die Temperatur 7 der Grundfläche z = Oals Ortsfunktion T(x, y, 0) 
vor, während für die Deckfläche T(x, y, h) = T, = konst. gelten soll. Die zugehörige Kon- 
vektionsströmung bestimmt sich dann durch das Differentialgleichungssystem (2.13), (2.15), (2.16) 
bei geeigneten Randbedingungen nachdem man in (2.15) die Lösung der Wärmeleitungsgleichung 
für dieses Problem eingetragen hat. In dieser Form dürfte esim Augenblick wohl fast hoffnungslos 
erscheinen, nach einer expliziten Lösung für w zu suchen, wenn man sich vergegenwärtigt, wie 
kompliziert allein schon der Aufbau der Lösung der Wärmeleitungsgleichung T7(«, y, z) sein 
kann. Wir sind also zu weiteren einschränkenden Annahmen genötigt, um konkrete Aussagen 
über die einsetzende Strömung machen zu können. Es ist klar, daß sich diese Annahmen auf die 
zunächst willkürliche Bodentemperaturverteilung beziehen müssen. Im Fall der zweidimensiona- 
len Strömung (x, z) wollen wir voraussetzen, daß die Bodentemperatur stückweise konstant ist 


7, 2 6 
ee 


Na.0), = | 


während für die Deckflächentemperatur nach wie vor 


tra) ul, 
gelten soll. | 
T, wählen wir hierbei so, daß die Temperaturdifferenz T,— T, dividiert durch h, —P, = (T,— Ty)/h 
größer oder gleich dem kritischen Gradienten der RAYLEIGH-JEFFREYSschen Theorie (= — Pre) 
ist. T, werde dagegen so gewählt, daß —ß, = (T, — T,)/h kleiner oder gleich dem kritischen 
Gradienten wird. Hierbei wird später die Differenz 7, — T, beliebig klein angenommen werden, 
jedoch so, daß die folgenden Ungleichungen erfüllt sind: 


— Pi 2 — Pirit » — PB, S — Panik. 
Bei genügend kleiner Differenz T,— T,in unserem Fall ist es dann naheliegend, mit folgenden 
Lösungen der Wärmeleitungsgleichung in (2.15) zu arbeiten: 
| T+M% .|e2]<b 
Ix;2) = 5 Pı el 
Tr Ba 20 1ctm.b,. 

Wir haben damit unser Problem auf den Fall eines konstanten — jedoch in den beiden 
Bereichen |e|<b und |x| > b verschiedenen — Temperaturgradienten zurückgeführt. Die 
Lösung bereitet jetzt keine Schwierigkeiten mehr. Selbst bei dieser starken Vereinfachung des 
anfangs gestellten Problems ergeben sich interessante Ergebnisse. In den Bereichen |x| <b und 


[x| >b hat man nun zwei grundsätzlich verschiedene Strömungen, die bie —= + b stetig an- 
einander schließen. In |x| < b liegt eine in x periodische Zellularströmung vom BENARDschen 


a 
Typ vor. Die Zahl der Zellen ist hier abhängig vom Verhältnis g- In |x| > b klingt dagegen w 


exponentiell mit — |x| ab. Hier liegen keine Zellen vor. Es handelt sich hier sozusagen um eine 
Sekundärströmung zu den Zellen, die bis zum Unendlichen reicht (siehe Bilder 2 und 3). Dieses 
Zusammenwirken zweier grundsätzlich verschiedener Strömungsfelder (periodisch und exponen- 


tiell abklingend) wird im Experiment häufig beobachtet und ergibt sich hier recht zwangsläufig. 
s+ 


u F eZ EIwer Er eur 


- u Be TEE, 
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3. Die Differentialgleichungen des Problems 


Das zähe, wärmeleitende Medium befinde sich über der x, y-Ebene. ‚Die Schichtdicke sei h. 
Die Beschreibung der Strömung erfolgt durch die 6 Variablen: Geschwindigkeit im = (u, v, w), 
Druck p, Dichte og und Temperatur d. Für diese 6 Größen bestehen 6 Gleichungen: 3 NAVIER- 
Storzssche Gleichungen, die Wärmetransportgleichung, die Kontinuitätsgleichung und die 
Wärmeausdehnungsgleichung (Zustandsgl.) 


Mut tn +rAu es Pe 0), 
Ruten 4-4 nit A (2.2), 
Mut tn = g+uAm 2 (2.3), 
ru t+ı +0 —kad a er ee He (2.4), 
+40 en (2,5), 
eier en (2.6). 


Die auftretenden Konstanten!) bedeuten hierin: » = kinemat. Zähigkeit, k = Temperaturleit- 
fähigkeit, « = Ausdehnungskoeffizient. 


Wir nehmen nun an, faß es sich bei der Konvektionsströmung um eine Störung des be- 
wegungslosen, stationären Wärmeleitungsvorganges mit der Temperaturverteilung T = T(«, y, z) 
handelt. Dementsprechend machen wir die Ansätze: 


9» —9,= Tl, y,2) + 0& 2, PEDPESE, ID EPEHZDTS FI 


wobei AT=0ist. © ist hierin die konvektive Zusatztemperatur und P der konvektive Zusatz- 
druck, und es gilt 


: | 
= —ya(ll—aT) a N 


Diese Ansätze tragen wir in (2.1) bis (2.4) ein und beachten, daß sowohl u, v, walsauch @ und P 
von erster Ordnung klein sind. Eine Linearisierung führt zu dem folgenden System: 


DIE 
ne 5 en. > 
ee er 
ee: +9 Ar : ti ara ae DR 
dw 1 oP 
7 De +9 09 F ya er ee Or 
Ro) oT oT oT 
rn =kA9—u Em Pr in (2.12), 
cu ov dw 
(= ii — 
ar + Em n EEE (2.13). 


Hierbei wurde die Dichte o als konstant vorausgesetzt, ausgenommen dort, wo sie mit g 
multipliziert wird und damit zu dem Auftriebsterm 9& O führt. Weiterhin ist zur Ableitung der 
Differentialgleichungen, wie man unmittelbar bestätigt, noch die Annahme 


|» To, 4 he en. ee 


erforderlich. Diese auch in der Rayveısuschen Theorie auftretende Bedingung besagt, daß das 
Produkt von x und d— 9, (im vorkonvektiven Zustand) viel kleiner als 1 sein muß. Im all- 


gemeinen ist diese Bedingung von selbst erfüllt, da & sehr klein ist (z. B. Paraffinöl 20° C,& = 6,7 
x. 10% grad}). ß i 


eh, ’ x 2 en i 
ER ER x &: f del. hier », k, « als konstant voraus. Bzgl. des Einflusses einer Variabilität von v verweisen 


ve 


a A 


np 
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Durch Elimination von P und © aus dem System (2.9) bi N 
Re y (2.9) bis (2.13) gelangt man zu der folgenden 


ee 0 oT oT oT 


tr (15): 
Hierin bedeutet: 

0° 0° 

20 dp 

Zu (2.15) treten noch’zwei weitere Gleichungen für u, v, w hinzu, damit die Bestimmung der Ge- 

schwindigkeit möglich wird. Durch kreuzweise Differentiation von (2.9) und (2.10) erhalten wir 


und P} e 


7) ou ov 
ra) (a 2) =: 0 ER rs Re are. (2.16). 
Wir nehmen dann noch die Kontinuität 
ou oD ow 
LE ee EC EN on (2.13) 


(2.13), (2.15) und (2.16) reichen zur Bestimmung von u, v, w aus. Es treten zu diesen Differen- 
tialgleichungen noch die Randbedingungen hinzu, die für eine horizontale feste Oberfläche 
folgende Form haben: 


wer m 0% REN BEER. (2.17). 
Für eine horizontale freie Oberfläche gilt dagegen 
ou OD 
im = 0% Ge Re (2.18). 


Wir merken an, daß die Bedingung © = 0 sich auch in w allein schreiben läßt. Auf Grund 
der oben durchgeführten Elimination ergibt sich die Beziehung: 


(4-4) 40=9540. 


ot 
Da © = 0 für alle x, y der festen oder freien Oberfläche ist, folgt 4,0 = 0 und damit 
0) 
en = Z) Eee (2.19). 


Damit bleibt das System (2.13), (2.15), (2.16) bei den oben angegebenen Randbedingungen 
zu lösen. Der Druck P und die Temperatur © werden dann anschließend bestimmt. 


3. Die Temperaturverteilung des vorkonvektiven Zustandes 


Betrachtet man das System (2.13), (2.15), (2.16), so erkennt man, daß hier entscheidend die 
Temperaturverteilung T = T(x, y, z) des vorkonvektiven Zustandes eingeht. Diese muß also 
zunächst als Lösung der stationären Wärmeleitungsgleichung AT = 0 bestimmt werden. Als 
Randbedingungen kommen hier etwa in Frage 


T&,y,0)=f(, y), RT Te (3.1a, b). 


Man sieht nun sofort, daß unser Konvektionsproblem nahezu hoffnungslos erscheint — 
insbesondere die Suche nach einer expliziten Lösung für w — wenn T(t, y, z) recht kompliziert 
wird. Wir wollen daher die Vorgabe (3.1a) möglichst einfach wählen, so daß einerseits das 
Wesentliche unseres Problemes erhalten bleibt, andererseits aber Aussicht auf eine explizite 
Lösung besteht. Wir nehmen daher die Bodentemperatur etwa als stückweis konstant an: 

Boelanelr = 
X, U, =], = . 
% Tr, al >b 


Ist nun die Temperaturdifferenz |T, — T,| genügend klein, so kann man bis auf Fehler der 
Ordnung O (|T,— T3,|) die Temperatur des Grundfeldes durch 


ae he 2|<b 
I2.9,2) = T n undalley .. . (3.9) 
Dee re mr Bz, lb 


h 
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ersetzen. Das geht auch daraus hervor, daß die Lösungen der Gleichung AT = 0 stetig von den 
Randwerten abhängen. Wir wollen diese Approximation des Grundfeldes nun für die Lösung der 
Differentialgleichung (2.15) benutzen. Hierbei ist jedoch zu beachten, daß in (2.15) der Gradient 
von T auftritt. Für |x] < b und |e| > b wird man (3.3) als Näherung für 7(z, y, z) auch 
noch in (2.15) nehmen dürfen. Schwierigkeiten bereitet allein die unmittelbare Umgebung 
von |z| = b. Wir gehen daher so vor, daß wir (2.15) in der vereinfachten Form (4.1) für 
Ic] < b und |x| > blösen, die Knüpfung der beiden Lösungen bei |c}—=b unter der folgenden 
Bedingung vornehmen’): 


HN 
a) i — = rn ° = — Brit +0, 


T—T, 
b) u ah a lea ine 


4. Die Lösung des Konvektionsproblems im zweidimensionalen Fall 


Nach den Ergebnissen des 3. Abschnittes vereinfacht sich unser Differentialgleichungs- 
system (2.13), (2.15), (2.16) erheblich. Für (2.15) erhalten wir jetzt 
> ‚ 


(5-4) (5 #4) Aw + 98 BıaA = 0 12) neu ER al), 


wobei hier und im Folgenden der Index 1 für |x| < 5 und der Index 2 für |x| > b eintritt. Es 
entsteht also eine Differentialgleichung für w allein, mit deren Lösung im stationären Fall wir uns 
jetzt beschäftigen wollen. Wir haben also: 


vykAAAw+goßaAw=0 ..:-........ (4.2). 
Hier machen wir den Produktansatz | 


w.= Fi(2) cf a rn RR IE ER (4.3). 
Die Variablen x, zlassen sich trennen, wenn wir verlangen 
u C 
May wi hs@. u... Se (4.4) 
mit C,, als Separationskonstanten. Für F(z) kommt dann 
2 Cara G 
vklaa a x Fo+go a EFQ-0 2.2.2... (8). 


Um den Sachverhalt möglichst einfach darzustellen, beschränken wir uns weiterhin auf den Fall 
zweier freier Oberflächen beiz= 0 und z=h. Wir betonen ausdrücklich, daß die im fol- 
genden zu entwickelnde Methode sich grundsätzlich auch beim Vorliegen anderer Randbedin- 
gungen durchführen läßt, wenngleich sie dann auch — wie schon in dem von RAyreıcH [2], 
JEFFREYS [3] und später von PELLEw und Soutuweız [6] behandelten Fall — zu mühsamen 
Rechnungen führt. Bei zwei freien Oberflächen erhält man als Lösung von (4.5) sofort F(z) = 


Erz ee } 
sin ==. Man bestätigt unmittelbar, daß F(z) die Randbedingungen (2.18), (2.19) erfüllt. Für 


das jeweilige Parameterpaar Cy., ßı, erhalten wir in diesem Fall die folgende Eigenwert- 
gleichung: 


2 v k(Cıa Pr 702)? 1 4 & Pıa Cya hi = 0 A Se lt LE rd (4.6). 
Führt man hier die folgende Abkürzung als RAyueıcH-Zahl ein 
ER 2.2 N, 
R=—- een ei" 
so geht (4.6) über in 
GC. —n)3 
R we ( 1,2 
1, C, RER 


ar 2 a Relation zwischen C,, und R,, muß sowohl für |c| > b als auch 
RR ie ; sein, damit eine nichttriviale Lösung w existiert. Sie ist eine notwendige und 
ne e Bedingung für die Existenz einer Eigenlösung der Form (4.3), (4.4),(4.5). Man spricht 

Shalb auch häufig von einem zweiparametrigen Eigenwertproblem. In Bild 1 ist (4.8) ne 


*) Die Symbole -+ 0 bedeuten wie üblich den links- oder rechtsseitigen Grenzwert. 
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tragen. Die Kurve R= R(C) besitzt für C< 0 ein positives Minimum und zwar bei Cih= 
— 70/2, Rmin = 27 nt/4. Dieses Minimum und auch der ganze Ast dieser Kurve für C<0O 
liefert mit (4.4) periodische Lösungen f(x), d.h. dieser Kurvenast entspricht den Zellular- 
strömungen. Der andere Kurventeil mit C >0 (nur R > 0 ist bei uns sinnvoll) liefert mit (4.4) 
monotone Lösungen vom Typ der Exponentialfunktion. Wir benutzen nun die Forderung (3.4), 


worin — krit denjenigen Temperaturgradienten bezeichnet, der zu Rum = 27 n“/4 gehört. Im 
Falle (3.4a) erhalten wir mit (4.4) 


Ze Pe 
1(&) + Im 1) —) R den Ih) = A cos % y Es er re (4.9a), 
R 
1000 
800 
27 _4 
Amin“ zT 
600 
400 
200 
“ Ui 2 2 
Cnn="7 u “T 
-20 E00 10 Ser 
-200 | 
1 
—4L00 ı 
—-600 


Bild1. Graphische Darstellung der Eigenwertgleichung (4.8) 


EZ 


Nr 5 ö 2 i > , Fä I aa 

, E 8, En. ET 
En) = 0, ha) = B 

BA AR A aa 
Insgesamt erhalten wir damit als I g von (4.2), ı 
x = + b verlangen und wegen der Symmetrie in x nur den 2 


RZ 
w= A 


(z—b) 
= Asin BE oa A 


h  Yh 
Wir wollen sogleich an ee Darstellung eine Betrachtung über die Ableitung = bi = 


een 
-, 


knüpfen. Die Stetigkeit dieser Ableitung würde verlangen \ = 
AR TE _ a Ab .. zu. le 
—— (b,z) = — — sin 7 si —— (b, z) = — —— sin — 0085 ——, pr 
Fa Y2h "pn ir 93 h h_ Y2h Ye 
d.h. . ’ 


sb _sR_ ee Ir Br, k 
ee Cr 


Dies stellt eine Bestimmungsgleichung für b dar und besagt, daß im allgemeinen die ° 
w-Verteilung bei = b einen Knick hat. Dies überrascht nicht und ist eine natürliche Folge * 
unserer Approximation der Temperaturverteilung T in (3.3). — An (4.10b) läßt sich auch leicht 
ablesen, wie rasch w, in x-Richtung abklingt. Es ergibt sich für eine Abnahme der Amplitude 
auf die Hälfte des bei x = b vorhandenen Wertes 


Br 1 
ee u 
e 3 B 
d.h. 
4x 1 0,6931 0,6931 
er n2= aa or a a (4.12) 


Wir gehen nun zur Bestimmung von u, , über. Hierzu verwenden wir die Kontinuitäts- 
gleichung (2.13) mit v=0. Es wird also 


x ER h h V2h 
und damit 


1,2) = —MWAcos — sin —— Tor Ri (z). 


Da aus Symmetriegründen u,(0, z) = 0 sein muß, folgt 9,(2)=0. Genauso erhalten wir 


an 2. u SE SZ 
&x E hr SIR 
und damit 
_3&-d 
ug(X, 2) = A Es 3 + 9(z) . 


Ö h Y2 h 


Die Funktion 9,(z) bestimmt sich aus der Forderung der Stetigkeit der u-Funktionen beixr=bzu 


922) = — A cos 7 (Y2sin 2 *2 + zes). 


h Rh 9 Yaı 
Damit erhalten wir die Darstellungen: 
a nz. nz 
un =—J2A cos zy su on elta (4.13a), 
@—b) 
_AE 052 zb (. nz. ab 
U; 5 008 7 cos Bu h 1)— 2A cos , sin a (4.13b). 


EFT GE NEQEN FOREN YATNRTANFRER 


wwrn'% 
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3 ; ou. 5 : 
Die Ableitung = ist bei x = b im Gegensatz zu = stetig, wie man mit der Kontinuitäts- 


gleichung unmittelbar bestätigt, da = stetig ist. 
z 
x Dagegen verwundert zunächst, daß im allgemeinen u, für x — 0 nicht verschwindet. Es 
ommt 


j nz 
lim u(, z) = — A cos — 
> +0 h 


TL sch sch 
— en en 4.14). 
3 "Don eu a) a. 


Dieser Ausdruck wird nur dann für alle zmit 0 <z < h Null, wenn gilt 


zb 7 
1 ES er a 4.15). 
V2h V2ö a 
Diese Beziehung mit (4.11) kombiniert, ergibt 6? = — n?/2, was wir aber ausgeschlossen 


haben. Wir können also selbst bei geeigneter Wahl von b höchstens eine der Bedingungen (4.11) 
oder (4.15) erfüllen. Daß im allgemeinen u, für x — oo nicht verschwindet, ist jedoch durchaus 
nicht verwunderlich, dies liegt vielmehr am zweidimensionalen Charakter der Strömung. Ein 
Blick auf Bild 2 bestätigt das, wenn man den Massenfluß zwischen zwei Stromlinien für |x| > b 
im Endlichen und Unendlichen betrachtet. Im rotationssymmetrischen Fall— den wir anschlie- 
Bend behandeln — wird sich zeigen, daß unter denselben Voraussetzungen alle Geschwindigkeits- 
komponenten im Unendlichen verschwinden. — Die Bestimmung der Stromlinien erfolgt am 
einfachsten über die Stromfunktion. Ausgehend von den Darstellungen (4.10a, b) für w,, und 
(4.13a,b) für u,, erhält man mit der Kontinuitätsgleichung (2.13) im ebenen Fall nach kurzer 
Zwischenrechnung 


V2 RES ICE 
FE) ee oT, EEE ER (4.168), 
az) 
ae en sin #2 sin =? = 2 (. Er ze . (4.16b). 
7 h Ph %ö Yen 


Wir haben in Bild 2 Stromlinien für die folgenden Werte von b gezeichnet: b = Y2 h/2, 
3 V2 h/4, V2 h, 4 V2 h/3. (6 = 2). Die Bezifferung der Stromlinien ist dabei so gewählt, daß für 
|z]|<b Y, folgende Werte hat: 


2 


a Fr 
ea re 01,0 005,0 
A h(— V2r) 2 
Für |x| > b hat ', folgende Werte 
SUR BAM =U; #015; +0,87 ...5, E10. 
A h(— V2/x) 


wobei das obere Vorzeichen für die ersten zwei Strömungen von Bild 2 zu nehmen ist, das untere 
Vorzeichen für die letzten zwei. 


Das Ergebnis kann man folgendermaßen zusammenfassen: Für |v| < b und |z| > b liegen 
zwei verschiedenartige Strömungen vor, die bei x = + b stetig aneinander schließen. Innen 
handelt es sich um eine periodische Zellularströmung, während außen die Geschwindigkeits- 
komponenten bzgl. x einen monotonen Verlauf zeigen. Die Periodizität der Innenströmung 
tritt hierbei natürlich in Abhängigkeit von der Wahl von b nur mehr oder weniger stark in Er- 
scheinung. 


5. Die Lösung des Konvektionsproblems im rotationssymmetrischen Fall 


Wir können uns hier kurz fassen, da die Lösung in diesem Fall analog zu derjenigen des 
ebenen Falles gewonnen werden kann. Die Temperaturverteilung des Grundfeldes sei dieselbe 
wie in (3.3) nur mit der Fallunterscheidung r < b und r > b. An Stelle von (4.3) kommt jetzt 
bei zwei freien Oberflächen der folgende Produktansatz: 


 — sin = fi.(r) An. a Pe 


Damit wird bei Rotationssymmetrie 


&w 1 cw ew 
Ber an 
or? Tor 02 


. (5.2). 


(Ka + 1 


> 
[ 
ea 
» 
—ı 
I) 
[7 
a 
un 
. 
B 
m“ 


0 ZB MEN 08.010 7 RR eee FaoN Bee 


P 


Bild 2. Quantitative Darstellung der Stromlinien im ebenen Fall. (Die Bezifferung wird im Text erläutert) 


An Stelle von (4.4) kommt jetzt 


nr v 


v r ig C | | 
them, 2 N 


mit C,, als Separationskonstanten. Tragen wir nun (5.3) in (4.2) ein und beachten (5.2), so er- 
halten wir dieselbe Eigenwertgleichung wie im zweidimensionalen Fall, nämlich 


(Ca — ga Bin GO 


a an ll nic a nn 
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Wir können also bei unseren weiteren Betrachtungen Bild 1 auch hier verwenden. Stellen 


wir nun wieder die Bedingungen (4.9 a, b), so kommen für (3.3) folgende Lösungen: 


Fi 1 , gu # GUT 
i+-htgah=0, dh O=Ad & A RR (5.52). 
Die Lösung N, (n r/ V2 h) entfällt, da sie für r— 0 nicht beschränkt ist?). Ganz analog erhalten 
wir im zweiten Fall 

BT Ar [08 RENT 

k+-h—-ah=0, d.h. ke)=BiH® ( ö ri) BRFETS (5.5b). 

Die Lösung J,(iör/h) entfällt, da sie für r— oo nicht beschränkt ist. Die Funktion f,(r) 

hat beiree!lem B die Eigenschaft reell zu sein und für r — oo nach Null zu gehen. Damit erhalten 
wir, wenn wir noch Stetigkeit beir =b verlangen 


IT GEL 
w, —Asin "2, © ze TER te ge a (5.68), 
ı b 
le 
Ann up (10) HER ER (5.6b) 


— 7 
BD: WR; ß 
ia (16, 


Wir wollen nun die Radialgeschwindigkeit U(r,z) bestimmen. Es bestehen die folgenden Zu- 
sammenhänge 
u= U(r,z)cosp, v = U(r, z) sin p 


mitz=rcosgy, y=rsing. Hierdurch ist bereits (2.16) erfüllt, wie man leicht bestätigt. 
Die Kontinuität ergibt 


cu 1 ev ee U U, 
oz ' . r 
Mit 
nz 

U,.(r, 2) = 1,.(r) cos Br 
erhalten wir für g,(r) die Bestimmungsgleichung 

; 1 7 

91,2 ir ee J1,2 + fıe ; U, (5.8) 


mit den Funktionen f},(r) aus (5.5a, b). 
(5.8) hat die Lösung 


er D h 
= | TE 69) 


D, . sind hierbei Integrationskonstante. Nach Ausführung der Integration in (5.9) erhält man die 
Ergebnisse 


u) =— WA, e = a Et ee (5.102), 


Sc". Sa ED 
g.(r) = BZ? Hi (154) ae rare (5.10b). 


Die Konstante D, muß Null sein, da sonst g,(r) für r — O0 nicht beschränkt wäre. D, kann 
so bestimmt werden, daß ein stetiger Übergang von g,(r) zu 9,(r) beir = b gewährleistet ist. Ins- 
gesamt erhalten wir also 


Un, 2) = — VBA cos = 7 © nl ee er Fe (5.11a), 
nn 2 RE LET 1, Da re) re (5.11b) 
27 0) - Has Fr En 


8) Die Funktionen J„, Nn, H% bedeuten hier und im Folgenden bzw. die BEsseusche, die NEUMANNsche, 
die Hankesche Funktion der entsprechenden Ordnung. 
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mit 
en 
Be A Te nn ba ana (5.12) 
un (16-7) 
und nz b 
SE 
D=--Ab|RA + T anwl[is” 
; | 7 Ir ne ioH (195) et (Gl 


: ed 
} iHW ( 6) n) 
Damit ist das Problem auch im rotationssymmetrischen Fall gelöst. U und w sind stetig bei 
N bund beide Funktionen verschwinden für r — oo. Letzteres leitet man sofort aus bekannten 
Eigenschaften der Besserschen und Hankeuschen Funktionen her, (siehe z. B. [7], S. 132). Auch 
in diesem Fall kann man leicht die Stromfunktion bestimmen. Es kommt nach Benutzung von 
(3.6, b) und (5.11a, b) mit der Kontinuitätsgleichung (5.7) nach kurzer Zwischenrechnung: 


v2 . nz ar 
A) =— has hl n) EEE WEM Eiepsohdke) = RER BR (5.14a), 
y2 nz ab 7 "(ur r b 
Pr) =——hAsin—|5J | = - ranp (07) enp (167) 


(5.14b). 
In Bild 3 haben wir Stromlinien für die folgenden b-Werte gezeichnet: b = by/2, 3 by/4, Das 


4bo/3. b, ist hierbei durch die erste nichttriviale Nullstelle von J, (r bo/ y2 h) —= (0 gegeben. Die 
Bezifferung der Stromlinien ist analog zum ebenen Fall gewählt. Für r < b hat /, folgende Werte 


yY 
Be: 170,125°9 20,253. 05, 
ar) 
IT 
während für r > b gilt 
1 
2 % 01-015 0 E05. 
rt) 
IT 


Das obere Vorzeichen gilt für die ersten zwei Fälle von Bild 3, das untere Vorzeichen für die 
letzten zwei Fälle. | 

Das entstehende Strömungsfeld zeigt gewisse Ähnlichkeiten zum zweidimensionalen Fall. 
Für r < b liegen hier Ringzellen [4a] vor, während für r > b die Geschwindigkeitskomponenten 
wieder ein monotones Verhalten zeigen. Auch hier ist es wie im zweidimensionalen Fall so, daß 
die Periodizität der Innenströmung in Abhängigkeit von der Größe b nur mehr oder weniger 
stark in Erscheinung tritt. Bild 3 zeigt handgreiflich, wie es zur Bildung der Zellularströmung 
für r<b kommt. Eine solche Kombination der zwei verschiedenen Strömungsfelder für r <b 
und r > bist oft im Experiment beobachtet worden. Dieserhalb verweisen wir auf die schönen 
Abbildungen in der Arbeit von v. TIPPELSKIRCH |5 c]. 


Herrn Prof. Dr. H. GöRTLER danke ich herzlich für kritische Bemerkungen zu dieser Arbeit. 
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Elastisch-plastische Deformationen von Platten 


Von A. LANGENBACH 


Die Durchbiegung w(x, y) einer elastisch-plastischen Platte wird durch eine nichtlineare partielle Diffe- 


rentialgleichung vierter Ordnung beschrieben. Randwertprobleme für diese Gleichung werden auf ein 
Minimum-Problem zurückgeführt und mit funktionalanalytischen Methoden gelöst. Die betrachteten Rand- 
bedingungen entsprechen der eingespannten, aufgelegten und freien Platte. 


The sag w(x, y) of an elasto-plastic plate is determined by a non-linear partial differential equation of 
the 4th order. Boundary-value problems related to this equation are reduced to a minimum problem and 
solved by the methods of functional analysis. The boundary conditions here considered correspond to the 


clamped, supported, and free plate, respectively. 


IIporu6 w(x, y) ynmpyro-nasctnyeckoß IIIACTUHKH OIMChIBAeTCA HeJIMHeÜHBIM YPaBHEeHMEM 
yerBeproro mopAmka. HpaeBpie 3ayaym MIA 9TOTO ypaBHeHuUA IIPHBONATCH K HEKOTOPOoH 
IIpoÖeMe MUHUMyMA MH peialTcH MeTolamu (yHKIMOHAJABHOTO aHasmaa. PaccMoTpeHHbIe 
KPaeBble YCJIOBUA COOTBETCTBYIOT 3AIleMJIeHHOÜ, ONEPTOÜ U CBOOOAHOA IIIACTHHKAM. 


1. Die Gleichungen 


Wir betrachten eine Platte, deren Last senkrecht zur Mittelebene wirkt. Unter der Ein- 
wirkung der Last erhält die ursprünglich in der x y-Ebene gelegene Mittelebene die Durchbiegung 
w(z, y) . In dieser Arbeit benutzen wir die geometrischen Relationen der linearen Elastizitäts- 
theorie, so daß wir die Deformationen in der Form 


ae A Da BR ee 

a 0x2 ’ Tu aa’ 

„_Ö__,® a, m, 

va ee öy? ’ Ye Te au ; . ID 
ow ev fan) 

= == A R a en en 

a er >> 


annehmen. Hierin sind ü(z, y, 2), d(, y, 2), W(x, y,z) die Verschiebungen in einem beliebigen 


= hesh £ 
Punkt (&,19,2)e 2x ee = des mit der undeformierten Platte identifizierten Zylinders. 


Die Formeln (1.1) unterscheiden sich nur durch die in der Plastizitätstheorie übliche Inkompressi- 
bilitätsbedingung e=&+ & + & = O von den bekannten Relationen der linearen Plattentheorie. 
Das Plastizitätsgesetz formulieren wir nach Hrxory durch die endlichen Beziehungen 


,—0—=4 g(H?) & 5 Toy = 2 g(H?) Yıy» 


oy—o=4g(H?)e,, Tg; = 2 g9(H?) y,; , ee 
4 —0=49(H3) s,, Tys = 2 g(H?) yyz 

| 1 1 

ee = (z+0%y+o)» 3 5a 


wobei 


2 2 2 2 2 
IP =yoyHt Ya t Yet 3 &- rg + — a]. 2... (1.4) 


eine bekannte Tensorinvariante bezeichnet. Wir wollen an dieser Stelle nicht auf die Bedingungen 
eingehen, unter denen das Plastizitätsgesetz (1.2) gültig ist und verweisen auf die einschlägige 
Es an en übrigen aber weitgehend willkürlichen 

n gestat g vers 'nartiger nichtlinear elastischer Probleme 
die — direkt oder indirekt — eine bedeutende Rolle in der modernen Plastizitätstheorie spielen. 


Die Spannunge enü h : et & ; . 
ee an > Tyz genügen den Gleichgewichtsbedingungen, in denen wir die 
asse ‘ vernachlässigen. Diese Vernachlässigung ist übrigens praktisch bedeutungslos, da 


ee 572 
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man eventuell vorhandene Massenkräfte (Eigengewicht), die senkrecht zur Mittelebene wirken, 
durch eine zusätzliche Last berücksichtigen kann. 


00% OTay 004% 

Br oy ee 

OTay 00, Ty; Er 

= a ee Te RE (1.5). 
OTyz OTyz 00; e 

Ze 2 ae 


Neun Gleichungen (1.2), (1.5) bestimmen bei gegebenen Randbedingungen neun Funktio- 
Denon ty dh D, U, 
Man sieht sofort, daß der Ansatz (1.1) unsere Gleichungen im allgemeinen nicht erfüllt. Man 
greift daher auf diein der Plattentheorie übliche Mittelwertbildung zurück. Hierzu führen wir die 
Biegemomente 


hj2 n]2 
Gaz=l.orzdz, a RS A ee (1.6), 
i —hj2 —h]2 
das Torsionsmoment 
h/2 
RI | Tau Zul MEET vr ER SERIE ITSIREE FARe (1.7) 
und die Scherkräfte - 
n]2 n]2 
N ewlrt.,de, Nu] myald2 nn nA Are (1.8) 
—h/2 —h| 
ein. Mit den Beziehungen (1.1) finden wir für die Intensität der Winkelverzerrungen (1.4) den 
Ausdruck 
o?w\2 Ow\? uw Aw Aw n 
er Re ON 29ER 
M2=4z ( a + | a + | 5 a + m Re (1.9). 
Somit lassen sich mit der Vereinfachung (1.3) die Spannungen o,, 0,, T,, durch w(z, y) ausdrücken: 
ee Nor 
rn ale EEE 
= 80) (tz a)” 
ow 1 &w 
= — 89ER) En + a EEE (1.10), 
or 
=4# 2 
rer y 
und dementsprechend auch die Biegemomente 
n w Po 2 
u. [72 > DaeHh: 
G2 D g*[H*%w)] & 77 at 3 
" ee a De (1.11) 
> w 2 
IR *[ 772 Set va 
re} [A w)) er A 2 
und das Torsionsmoment 
D nn ow 
gt ZDNet  HE aee 1-12 } 
R-— ze] (1.12) 
wobei 
2 hl2 
N Eden ana ©. (1.13) 


—h/2 


an die Stelle der früher eingeführten Materialfunktion g(H?) tritt. Durch Integration über z 
findet man aus den Gleichgewichtsbedingungen (1.5) mit Rücksicht auf (1.10) 


2=h/2 z=h/2 / 
ee N te rin si se (1.14) 
und 
oN, ON, 
- Fe ARE ee 1.15), 
el 2y + pP y) (1.15) 
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wobei [en — p(x, y) die Gesamtbelastung bezeichnet. Wir multiplizieren nun die Glei- 
nn 


chungen (1.5) mit z und integrieren wiederum. Die ersten beiden Gleichungen liefern dann 


0G, oK 
a 
3K Pe ee (1.16). 
ee Si! 
Nena 
Setzt man (1.16) in (1.15) ein, so findet man die Momentengleichung 
0°G, OK 8G, 117 
Fr] ran öy? + p(&, y) ee a A. ( )» 
die in ausführlicher Schreibweise mit den Beziehungen (1.11), (1.12) die Form 
Ere. (dw 1&w a NER . 
rer 
0° o., (Ow 1 &w 1 
— \gHHY) I — + <=) PS U) -»:: 2.22.20. 1.18 
at | De (1.18) 


erhält [1]. 

Die Materialfunktion /(&) = 9*(&) E wird gewöhnlich experimentell ermittelt. Unseren 
Ausführungen liegt ein Diagramm zugrunde, welches für die Spannungsverfestigung charakte- 
ristisch ist (vgl. Bild). Wir weisen auf folgende Eigenschaften der Funktion I(£) hin: 


Für & > O gilt r 
a) E =D DE>R>0, re cnskh 


b) die Krümmung ist nicht positiv, 
c) g*(0) > 0. 
Aus den Eigenschaften a)—c) folgt 
d) g*0) Sg) > x 
denn S 
1: are, 
Eee aaa Seh 
Ne a] | 
wegen b) und 4 j 
Ne) 1 fdll) 
KERN —_ er 
IH) E E | Fra; 


ee nr DER 


IV 
X 
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Zur Lösung von Randwertproblemen für die Gleichung (1.18) benutzen wir ein allgemeines Ver- 
fahren der Variationsrechnung [2]. 


Der Gleichung (1.18) ordnen wir das Funktional 


a0) 
9 
Dlw) — | | vo de ppw U Se ee er 


Q 0 
zu. Wie die Gleichung (1.18) wurde auch das Funktional (1.19) zuerst von L.M. KATSCHANOW 
angegeben. Als Definitionsbereich Da sei die lineare Mannigfaltigkeit M derin @ + S viermal 
stetig differenzierbaren Funktionen erklärt. Weiterhin sei P&, y) € L,(2), wobei mit L,(2) der 
Raum der mit ihrem Quadrat in 2 integrierbaren reellen Funktionen bezeichnet wird. An- 
genommen, P(w) besitze ein Extremum auf einer linearen Teilmenge M, < M; dann gilt 
e DP(w, + en] r = 2 [tar H?%(w,, h) — ” p | I2=0,.,% (0) 


für jedes he M,, wenn wy(z, y) die Extremalfunktion bezeichnet. In (1.20) ist Aw, h) die bi- 


lineare Form, die der quadratischen Form H?(w) entspricht: 
Puh woh 1 Aw oh euch Aw Ah 


1 
Om? 0a | Or oe 2 0 OD ae ey Ox öy 
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H%w, h) = 
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‚ Zur Umformung der Variationsgleichung (1.20) benutzen wir die bekannten Relationen zwischen 


Spannungen und Deformationen am Rande der Platte: 
N=N,cosp+N,sing, 
G = G,008?p + G,sin®9 + Ksin2p, 


1 : 

HZ 5 (6) sin 2 +Kcos2p 
wobei @ den Winkel der äußeren Normale an den Rand S mit der positiven Richtung der 
x-Achse bezeichnet. Die Biegemomente G,, G, und das Torsionsmoment K stellen wir uns dabei 
nach den Formeln (1.11), (1.12) durch die Durchbiegung w(x, y) ausgedrückt vor, während man 
für die Scherkräfte N,, N, entsprechende Ausdrücke aus den Gleichungen (1.16) erhält: 


(7) m a 7/0Ww 1 m 0 n w 
Net te 


Neu D 5 n (rL&w)] N ey 5 [prLA®w)] ( a 3 a u 


Durch partielle Integration des Ausdrucks!) 
[g* H:(w, h) dQ2 
2 


erhält man 


jr H:(w, h) d2 = I% (w- + > a) h; cos“ + (or + D w..) h,cosny + W,yh, eosny] ds 


S 
[) 1 
en] 
S 


„Wrunloosnz]na+ |Punan. 


1 
cosntc+ = Io* (tor - 75; 2) cosny 


— n 
ö 

2 
Mit den Bezeichnungen 


0 [7 ö 
cosnx=coSsp, I N a N 
: 0 0 . ö 
cosny=S1NY, a re 


findet man 


a 


n A 1 El, . ah Eos 
jr H?:(w, h) do == g I R 7 u) | 25 sın (07) [Mer (07 _ Zn 2 
Ss 


1 oh; die 2 sin? 
+ (v1 -F Do, 4.) v3 sın 9 C0Ssp en 2 
+ Way = sino + "sin peosg ds 


d 1 
sinpcosp+ En " (ter 4 5 a) cos? 


GR 1 
-[Faleteta) 
5 


do 1 
fr ” (vv E 5 Wr) 


08 


in2 
sin? 


7) 1 
sinpcosp+ an " (vv En D :) 


e ° i 
Ar = [g* w,,] co®?p + An [9* w;,] sin @ cos | hds + I? whd2. 
s 


LP} 


\ i des Gaussschen Satzes zuläßt. Überdies 
1 hmen an, daß der Rand 8 von 2 die Anwendung des 
sei 2 ee sternförmiger Gebiete im Sinne von [3]. Bei Körpern, die man als Platten anspricht, 
dürften diese Voraussetzungen immer erfüllt sein. 


g 
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Nach einer Umstellung der Glieder ergibt sich hieraus mit den Bezeichnungen (1.22), (1.11), 
(1.12) und (1.23) 


x oh öh oH 
D [vr ii. nan = — [no [oma + are 
5 


8 8 
= ai °G nds+D | PwndQ 
+ [rar [rat [mna+ w 
j 8 5 8 2 b 
oöH _0G oh 
ae Sri: — — D|Pwhd2, 
-/b ER +, | Rd jea®+ | w 
5 Pr 


S 


da Integrale der Form j nn (H h) ds verschwinden. 
u. * * 
Ferner findet man leicht mit Hilfe von (1.11), (1.12) und (1.23), daß zwischen den Funktio- 


nen N, G und H die Relation ’ 
"ah 06 


besteht. Endgültig erhält damit die Variationsgleichung (1.20) die Form 


1% ” 1 
nn n) = | le Bw, — — n| d2 
2 
Sue h T a 1.24 
5 a + ra 5 [me + [Ps pP 2.1.24) 
Ss Ss [P} 


für jedes he M,. w(z, y) bezeichnet die Extremalfunktion, G, = G(w,) usw. 

Um das Verständnis der folgenden Ausführungen zu erleichtern, fassen wir die Ergebnisse des 
ersten Abschnitts noch einmal zusammen: Die elastisch-plastische Deformation einer Platte läßt 
sich durch die Wölbung der Mittelebene w(z, y) beschreiben, die einer nichtlinearen partiellen 
Differentialgleichung vierter Ordnung genügt (1.18). Diese Gleichung muß bei entsprechenden 
Randbedingungen integriert werden. Für den einfachsten Fall der eingespannten Platte erhalten 
wir:z.B. das folgende Randwertproblem: 


Pui= 5 P@, Yy) = nn an a a 
dw & 
uls=0, Se 0. Sr wa Dee 


Ss 

Eine geschlossene Lösung dieses Problems kann nur in äußerst einfachen Fällen angegeben 
werden. Daher gehen wir zum Minimum-Problem für das Funktional (1.19) über. Einerseits 
liefert dieses Verfahren später den Existenz- und Einzigkeitsbeweis für eine passend verallgemei- 
nerte Lösung des elastisch-plastischen Plattenproblems, andererseits bietet es praktisch unbe- 
grenzte Möglichkeiten zur Konstruktion von Näherungslösungen. 

Die Variationsgleichung (1.24) zeigt, daß eine hinreichend glatte Lösung des Minimum- 
Problems der Gleichung (1.18) genügt. Diese Schlußfolgerung wird im nächsten Abschnitt für 
einige Beispiele mit den in der Variationsrechnung üblichen Methoden begründet. Es besteht 
allerdings die Möglichkeit, daß das Minimum-Problem für das Funktional (1.19) eine Extremal- 
funktion w,(@, y) liefert, die keine vierten Ableitungen in @ + S besitzt, da diese in dem Funk- 
tional (1.19) überhaupt nicht vorkommen. Hierzu sei bemerkt, daß man die Theorie der elastisch- 
plastischen Platten wie auch die lineare Plattentheorie direkt auf Variationsprinzipien aufbauen 
und somit auf die Gleichung (1.18) gänzlich verzichten kann. Eben mit Hilfe eines solchen Prin- 
zips wurde das Funktional (1.19) von L. M. KATSCHANOw gefunden, während esin der Arbeit [4] 
durch rein funktionalanalytische Methoden aus der Gleichung (1.18) abgeleitet wurde. 

Das Problem (1.18), (1.25) wurde in der Arbeit [4] übrigens ausführlich behandelt. Dagegen 
konnten die sogenannten natürlichen Randbedingungen 


oH 
Zt N=0,. Gel 2 


die wir ebenfalls der Variationsgleichung (1.24) entnehmen, dort nicht untersucht werden, da sie 
wesentlich nichtlinear sind. 


in 
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Wir wenden uns nun Existenz- und Einzigkeitsfragen einer Lösung des Minimum-Problems 
für das Funktional (1.19) zu. Die Bestätigung, daß die von uns gefundene Extremalfunktion 
wirklich so glatt ist, daß sie die Gleichung (1.18) erfüllt, bleibt späteren Untersuchungen vorbe- 
halten. Gewiß ist nur, daß die Gleichung (1.18) bei entsprechenden Randbedingungen überhaupt 
keine Lösung besitzt, wenn unsere Extremalfunktion die hierzu notwendigen Eigenschaften nicht 
hat. Den Beweis dieser Behauptung findet man in [4]. 


2. Untere Schranken für das Funktional @(w) (1.19) 
Angenommen, das Funktional 


[a] 


H’(w) 


N / g*(E) dE dQ 
genüge auf M, der Ungleichung 
Pia a ne u (2.1) 
Q2 


für eine passend gewählte positive Konstante y. Dann besitzt ®(w) eine untere Schranke auf M,. 
Wegen der SohwaArzschen Ungleichung gilt nämlich 


Du) = Fa) — pl ll = el? — Ip Ill = — IR, 
lol =[fpraope. 


Aus der wichtigen Relation (2.1) folgt ferner 
2 


> a 
Du) = Fo) — 77; ll FW). 
Daher ist ®(w) positiv für 
‚„ #lpl? 
Fi) > ap —K 
Auf M, führen wir nun die Metrik des Raumes W@(2) [3] ein, die wir in der Form 
2, |ewle Aw|? &w ||? 
\wliwe — | ap ne en a (2.2) 


schreiben können, wenn wir die Ungleichung (2.1) und die Eigenschaft d) der Funktion /'(&) 
berücksichtigen. Wegen eben dieser Eigenschaft gilt auch 


F(w) > x | A&w) dQ =, wire 
2 
Es genügt daher, wenn wir das Minimum des Funktionals @(w) auf denjenigen Elementen von M, 
suchen, die in der Kugel 
RK 


Le ee ee (2.3) 


des Banacnraumes W@)(.2) liegen, da außerhalb dieser Kugel @(w) offensichtlich positiv ist. 
Wir betrachten einige Beispiele für lineare Mannigfaltigkeiten, M,< M, auf denen die 
Ungleichung (2.1) erfüllt ist. 
I) Wir erklären als M( die Mannigfaltigkeit derjenigen Elemente von M, die auf S ver- 
schwinden: 
DEREN RES OR Tr 


Für die Elemente der Mannigfaltigkeit M® gilt die Frieprıchssche Ungleichung 


ow\2 ow\2 1 2 IE 
&®) en er) dQ >= ” Ik RE ARE IERTEEE ER PN 
2 


2 
RAT. cw cw 
in der a die Seite des Quadrates bezeichnet, das Qeinschließt. Auf —, Fn wenden wir nun die 


Poıscar&sche Ungleichung an: rn 


dw\2 o?w\2 ow \2 "s Ow 2 

a ee B un 
er a2= 4) +) |? (} = 
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LP} 


9* 


E 2 Br n- nice FE Br + 
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Wegen (2.4) gilt aber 
[= [eoro=0. 


[0 5 - 
cw h i 
Daher erhalten wir mit einer entsprechenden Abschätzung für Pr] aus (2.5) die Ungleichung 


1 i d2w\? Aw \? (wi? 1Q 
cf" ins | + + ee] 
Q2 Q 
mit einer bestimmten positiven Konstante C. Die Ungleichung (2.1) folgt nun wieder aus der 
Eigenschaft d) der Funktion I'(&): 


C 


2 


Auf MW besitzt D(w) also eine untere Schranke. Angenommen, ®(w) nehme sein Minimum auf 
M® wirklich an. Dann gilt also die Variationsgleichung (1.20) für ein w) € M{) und beliebige 
heM@). Bezeichnet man mit M® diejenige Teilmenge von M{P, deren Elemente auch noch 


der Bedingung 


rw) 2, [wa0. s 


'h 

Re 
on |, 

genügen, so gilt die Gleichung (1.20) und mit ihr (1.24) erst recht für ke M®. Mist dicht in 

L,(2) und daher 


1 
Pw=zPp- 


Die üblichen Überlegungen der Variationsrechnung liefern außerdem die Randbedingung 
GuV) = 0, 
während die Zugehörigkeit von w(® zu M® die Randbedingung 
ul; = 0 
erfordert. w(® beschreibt also die elastisch-plastische Deformation der aufgelegten Platte. 


II) Das Minimum-Problem für ®(w) kann auch auf der linearen Teilmenge M® von MW 
betrachtet werden, da auch hier natürlich die Ungleichung (2.1) gilt. Nimmt man an, daß w(@) e M3 
das Minimum von ®(w) auf M® realisiert, so findet man 


1 PRO) 
Puw® = D han w@ls == Ur än = 
Ss 


w®) beschreibt also die elastisch-plastische Deformation der eingespannten Platte. 
III) Wir erklären nun M® als Teilmenge derjenigen Elemente von M, die den Bedingungen 


fwaQ=d; [en dR WW; Tue Ti (2.6) 
[p} 


Q Q 
re 


erklärt. Aus der Beschränktheit des Einbettungsoperators [3] WPXQ2)— L,(2) folgt wegen (2.6) 


1 22 e a 
— [wdQ2 < Fw Prw rw 
y ay? 0x &y 

{pP} 


genügen. Auf M® ist die Norm des W®(2)-Raumes in der Form 


n 2 2 2) 1/2 
\wllwen, ın = | | van -+ BE » a0 -+ | [es al + | 
[4 S S 


P] f] 2 


12 


ww 


Fw\?  ||w 
0x 0y 


|” 


2 2 


+ +2 


Br 


Die Ungleichung (2.1) ist also wieder eine Folge der Eigenschaft d) der Funktion IE). 
Fordert man, daß auch die Plattenbelastung Bedingungen der Art (2.6) erfüllt, 


[pa2 =); BU0 0% {py dB =Rin2E ge 2.7); 
2 £P} 
so folgt aus der Annahme, daß w(® das Minimum des Funktionals D(w) auf M® realisiert, 
1 Er} 2H(w®) 
Pu D Ds G(wP)|s =(; + N(w®)|s =0. 5:28 
s 


En, 9 
RN 
r 
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Die Bedingungen (2.8) bedeuten, daß am Plattenrand keinerlei Kräfte wirken. w®) be- 
schreibt also die elastisch-plastische Deformation der freien Platte2). In diesem Fall darf natür- 
lich die Belastung weder eine resultierende Kraft noch ein resultierendes Moment erzeugen, wenn 
sich die Platte im Gleichgewicht befinden soll. Daher sind die Bedingungen (2.7) auch notwendig. 
w%) erfüllt die Bedingungen (2.6). Für die Lösung der Gleichung Pw= p/D mit den Rand- 
bedingungen (2.8) ist das belanglos, da ıw(x, y) ohnehin nur bis auf eine lineare Funktion in x und y 
bestimmt ist. 


3. Die Lösung des Minimum-Problems 


D(w) besitze auf der linearen Mannigfaltigkeit M, eine untere Schranke; d sei seine 
untere Grenze. Überdies genüge ®(w) der Bedingung (2.1). Wir betrachten eine Minimalfolge 
{w„} < M;: 

Um, Don) di) EN RT (3.1). 


N— 00 


Ohne Einschränkung der Allgemeinheit können wir annehmen, daß die w, der Kugel (2.3) an- 


gehören. Wie in [2] bilden wir das Funktional 


1 1 w w 
000 1) = z Pl) +5 Du) — D 
Wegen (3.1) gilt 
0o(w,, w,) er Or(w,, 1.) >e0R 
pP, > oo 
Nach einfachen Umformungen erhält man 
14 
1 ‚[r s+t 
Or(wy, w,) >= 2 I} B K (® + D- (wy 2 wo) 
200 
da t dz s+t 
u iz: (2 - = (w, u) A H? (1 E- ZT (u u) 
um t d(dz t 
gr % [1 r en wu) mr e: Me (u. P2. = (w, — o)))| ds dtdQ. 
Führt man die Bezeichnungen 
Ss+4 
w + — = uw*, w—ıw,=h 


ein, so bemerkt man unschwer die Beziehungen 


d > des a 
— H%w*) = — H%“w*) = H%(w*, h 
Pur) = — Hw*) = Hw*, h) 
und ; 
did N m 
Ban! ey 
ae) N) 


Benutzt man die Cavoaysche Ungleichung in der Form 
[A&w*, n)® < H%xw*) H:(h) 
und die Ungleichung g9*’(£2) < 0 (Eigenschaft b) der Funktion J(£)), so findet man die Abschätzung 


W — W, 


13% 
Orly, w) > | | \ {2 9° [we] Axw*) + gr[A&wr]} R> Be ds dt dQ 
200 
und endgültig 


ort up =» |, 


2 


a Be Er (3.2) 


wegen der Eigenschaft a) der Funktion J'(£). 


2) Man bemerke, daß w{?? wegen der Bedingungen (2.7) das Minimum von ®(w) auch auf M realisiert. 


Kleine Mitteilungen - ar 
EEE ERBE EEE EBRERGRE se sn nn, 
i i ie Mini in der Metrik W@(2) konver- 
Die Ungleichung (3.2) bedeutet, daß die Minimalfolge {w,} in w; 
giert. Ihr ek sei ,. Offensichtlich liegt w„in der Kugel (2.3). Aus der Eigenschaft d) der 
Funktion I(£) folgt die gleichmäßige Stetigkeit des Funktionals F(w) und damit auch von P(w) 
in der Kugel (2.3). Daher ist 
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lim dw) =dw)=d .:... ch 


n—00 
i i i is fü Ö betrachteten 
Die Relation (3.3) kann man als Existenznachweis für schwache Lösungen der t 
Probleme der Plattentheorie auffassen. Offensichtlich ist ıw, eindeutig bestimmt; hiervon über- 
zeugt man sich durch Vereinigung verschiedener Minimalfolgen. 


Wesentlich für die vorstehenden Überlegungen sind die Eigenschaften a) und d) der Funk- 
tion I{&). Die Eigenschaft b) wurde nur für die Abschätzung (3.2) gebraucht, die auch mit den 
Eigenschaften a) und d) zu beweisen ist. Andersartige Bedingungen wurden in der Arbeit [4] 
untersucht. Funktionale der Form (1.19) treten häufig in der’Theorie der elastisch-plastischen 
Deformationen auf, so daß die Ergebnisse dieser Arbeit auch über die Plattentheorie hinaus An- 


wendung finden. 
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KLEINE MITTEILUNGEN 


Kriterien für die exponentielle Stabilität der 
Nullösung eines Systems linearer Differential- 
gleichungen mit veränderlichen Koeffizienten 


In einer Arbeit von CHARASACHAL [3] wurde auf 
ein System von zwei linearen Differentialgleichungen 
mit veränderlichen Koeffizienten ein Verfahren an- 
gewendet, das gestattet, hinreichende Bedingungen 
für asymptotische Stabilität der Nullösung aufzu- 
stellen. Die Methode für die Konstruktion einer 
Lsarunowschen Funktion in dem gegebenen Fall ist 
ursprünglich von CETAEv, aber auch von Razu- 
MICHIN und LEBEDEV und vielen anderen Autoren 
hauptsächlich zur Ableitung hinreichender Bedingun- 
gen für asymptotische Stabilität benutzt worden, siehe 
z.B. [2], [5], [6], [7]. In der vorliegenden Mitteilung 
werden die von ÜHARASACHAL aufgestellten Be- 
dingungen korrigiert [8] und auf den Fall von n 
Differentialgleichungen verallgemeinert; aus bekann- 
ten Sätzen wird dann abgeleitet, daß sogar exponen- 
tielle Stabilität besteht. Zuletzt wird an einem Bei- 
spiel gezeigt, daß diese Bedingungen für exponentielle 
Stabilität nicht notwendig sind. 

Es sei das System linearer Differentialgleichungen 
gegeben: 
de; E 
ZT Mer t tanz; i=1,2,...,n(l), 
in Matrizenschreibweise: 


dx 


ir Re ) 


mit A= (ar), = fa, %,.2.,%n). Die Ele- 
mente a;; der Matrix A seien samt ihren ersten Ab- 
leitungen in dem Intervall t=t, Z0 beschränkte 
Funktionen der unabhängigen Veränderlichen t: 
air = Airft). 


Wir betrachten die Matrizengleichung 
ATB+BA=O SE. Rule 


Die Matrix B sei symmetrisch: B= BT (BT ist die zu 
B transponierte Matrix); dann ist auch C notwendig 
symmetrisch. Auch das Umgekehrte gilt, wenn (3) 
eindeutig lösbar ist; das wird im folgenden voraus- 
gesetzt. Bilden wir aus den Elementen b;z = br; und 
Cik = Cki der Matrizen B und (© die Vektoren 


D8) = {b,1> Drss Dre b, ns Ds, Ds .. ., Dan} ’ 


Os — {C» Gas +» > Can} D 


so ergeben sich aus (3) n(n + 1)/2 lineare Gleichungen 
für die ebensovielen Elemente b;x. Diese Gleichungen 
schreiben wir etwa in der Form 


AB: 


Hat die Koeffizientenmatrix A die charakteristi- 
schen Zahlen A,, Ag, » . » „An (die stetige Funktionen von 
t sind), so hat die Matrix A bekanntlich [4] die cha- 
rakteristischen Zahlen ; + (L<SiskS<Sn). Da- 
her ist die Determinante 


F ‚den Ku 


., Cın Ca9» Cay » ..,Cons« 


di=' IE aA) 5 een 


1sisksn 


4A ist also regulär, wenn keine der Summen A; + Ak 
verschwindet (insbesondere also, wenn jede Zahl A; . 
einen negativen Realteil hat). In diesem Fall kann ° 
man in (3) die Matrix (© beliebig wählen, und die 
Matrix B ist dadurch eindeutig bestimmt; denn (4) 
ist, eindeutig lösbar. Nach der Cramerschen Regel 
läßt sich das Ergebnis in der Form 


_ Ak, 


bir dl ’ 


lSi<sk<n-::.(6) 


schreiben. 
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Die Determinante |A| steht in engem Zusammen- 
hang mit der Hwrwırzschen Determinante H,„ der 
Gleichung |A E — A| = 0. Das sieht man, wenn man 
die Formel von ORLANDo [4] - 


n(n+1) 
2 u At)... Mm 


1<i<k 


HR=(—1 


heranzieht. Demnach gilt wegen (5) 
n(n-+1) 
(1) ° |4|= 2% Hn ES), 
Die Elemente a; = a;;(t) von A sollen nun so be- 
schaffen sein, daß. für jeden festen Wert > i, die 
charakteristischen Zahlen der Matrix A negative Real- 
teile besitzen. Dafür ist notwendig und hinreichend, 
daß alle Hurwırzschen Determinanten der Glei- 
chung |A#— A| = 0 zu dieser Zeit t positiv sind. 
Diese Bedingungen schreiben wir in der Form 


5204,20, 44H Fürs ei). 
Dabei gilt 
ES ANFR 05 0 ) 
De. ) 
Net 1 = Re ) (10), 
) N) N er N 


wo % Ar die Summe der Hauptminoren k-ter Ordnung 
der Matrix A ist. Man weiß, daß in diesem Fall die 
Gleichung (3) jeder negativ definiten Form x7 C x eine 
bei festem t > t, positiv definite Form zT B x zuordnet. 
Wir werden unten in (13) speziell =—2H an- 
nehmen. 
Wir nehmen nun als Lsaruwowsche Funktion die 
Form 
n(n+1) 
v»v=(—1) ? |4laTBe al) 


und bilden ihre totale zeitliche Ableitung für die Diffe- 
“ rentialgleichung (2): 


d n(n-+1) d 
un Fe 2 
2) [san ame 02) 
wobei B die Gleichung 

MEBEBAZ—9IE 13) 
erfüllt. 


Nach dem Satz von SYLvESTER ist die Form (12) 
bei jedem ? > f, genau dann negativ definit, wenn die 
Hauptabschnittsdeterminante M7, der Matrix 

n(n+1) d 
Mi) aaa] 00 


für > ti, positiv sind: 

De NE Oral, > 0 für 1 >n,.(1d). 

. Die Elemente der Matrix M haben wegen (6) die Form 
nur | 

me=(—1) ? & |Al dir Ai) (16), 


wor Mes Ger 1 für ik 9 = Äfür VE% 
ist. Im übrigen ist die quadratische Form v gleich- 
mäßig klein, da die Koeffizienten, die dort auftreten, 
als Summen von Produkten der a;;(t) für =>, be- 
schränkt sind. 

Die Bedingungen (9) bzw. (15) sind nicht, wie es 
von: CHARASACHAL angenommen wurde, hinreichend 
dafür, daß die Form v in (11) positiv .definit und die 
= in (12) negativ definit ausfällt [8]. Das 
zeigt schon das einfache Beispiel 


yet y 


Form 


w—erig, Ein): 


- 
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Die Differentialgleichungen (17) besitzen das Funda- 
mentalsystem von Lösungen: 


@ er ) ; [0 il : 
e e 


Die Nullösung ist daher nicht asymptotisch stabil, 
obwohl, wie man. leicht sieht, für (17) sowohl die 
Bedingungen (9) als auch die Bedingungen (15) gelten. 
In diesem Fallsind die Formen v und Sr nicht positiv 
definit bzw. negativ definit. 


Wir sorgen nun dafür, daß die Form » positiv 


definit und die Form = negativ definit werden, 


indem wir an Stelle der Bedingungen (9) und (15) die- 
folgenden Bedingungen verlangen [1]: 


EN Ten Moschee nel, ee; 
fürt>t, wre), 
mM, >.05,.2195 2105 Mae 020, 
fürt>t, . (19), 


wo ö eine beliebig kleine, aber feste Zahl ist. 


Wir zeigen, daß die Nullösung von (1) exponentiell 
stabil wird, wenn die Koeffizienten a;z(t) samt ihren 
Ableitungen ersten in Zukunft beschränkt und die Be- 
dingungen (18), (19) erfüllt sind. Dazu zitieren wir 
zwei Sätze aus [5]. 


Satz I (von Pexsıpskr). Wenn für eine lineare 
Differentialgleichung 


2 Ma al) 


eine gleichmäßig kleine, positiv definite Funktion v 
existiert, deren für (2) gebildete Ableitung © negativ 
definit ist, so ist die Nullösung exponentiell stabil. 


Satz II (von MAtkın). Wenn die Nullösung von (2) 
exponentiell stabil ist, kann man eine gleichmäßig 
kleine, positiv definite Funktion v(z,t) so konstru- 
ieren, daß ihre Ableitung © gleich einer willkürlich 
vorgegebenen negativ definiten Form — w(x,t) der 
ganzzahligen Ordnung m mit stetigen und beschränk- 
ten Koeffizienten ist, und zwar ist 


v(z, t) = Sw(p(T,%, t), T) dr (20). 
i 


Dabei ist p(7, ,t) die im Zeitpunkt # vom Punkt x 
ausgehende Lösung. (Die so konstruierte Funktion v 
ist offenbar auch eine Form m-ter Ordnung.) 


Unter den Bedingungen (18) und (19) genügt die 
Form v allen Voraussetzungen des Satzes I; daher ist 
die Nullösung exponentiell stabil. 


Daß die oben ausgestellten Bedingungen nicht für 
exponentielle Stabilität notwendig sind, zeigt folgendes 
einfache Beispiel. Die Differentialgleichung 


EN 


hat die allgemeine Lösung 


= 2 eg ttasint r (22). 
Für (22) besteht die Abschätzung 
Ix| < |z,| el"! e=! aaa) 


Daher ist die Nullösung von (21) exponentiell stabil 
für jedes «. Aber für |@| > 1 ist die Bedingung (18) 
nicht erfüllt. Für (21) kann man nach Satz 1I die 
folgende Form als Lsarunowsche Funktion bilden: 


co 
v—_2r eure) dt... (24): 
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* Man sieht leicht, daß 


[00] 


1 <2 et ssint) e2(esin?—r) dt < ealel 
ealel 2 


D 
ist. Daher ist v positiv definit und gleichmäßig klein. 


Nun ist = — — 2x2. Die Form (24) genügt also allen 
Voraussetzungen des Lsarunowschen Satzes über 


asymptotische Stabilität. 


Wir berechnen die erwähnten Bedingungen für die 
Differentialgleichung 


dx der da E 
Te + (25). 
Das äquivalente System der Form (1) lautet 
&=y; y=2; 2=—r—al)y—z (26). 


Die Hurwırzsche Determinante H,„ für die cha- 
rakteristische Gleichung von (26) ist 


a) 
7991| =g—1% 
07021 
Die Bedingungen (18) reduzieren sich auf 
at) — EI Art. 32 (27): 


Die quadratische Form (11) nimmt die Form an: 
v=-ß-—I+MFH+AarI LM Yr@g+DdIZ 
+22g +l)ay+2(—1)x2+6ya. 
Dann ist 
= [g(29—1)—2(9—1)]2?+ [9 (29 +1)—2(g—1)]y? 
+B9—2@— 1] +4g92y+2gxz, 


wobei qg = q(t),q 4 gu) ist. 

Die Bedingungen (19) lauten 

M=24—-1)—-429—1)>0, 

M, = 2 —)— 429g —1)][2(g—1)—4(2g+1)] 
— 4,0, 

1,=-[2@—-1)—-91;—Pl2@—h)—42g-+1)] 
0: 


Diese Ungleichungen sind bestimmt dann erfüllt, 
wenn 


ee Buben 1 


wobei &(ö’) > 0 ist und &(ö’) — 0 für Ö° — 0. 
Es sei nun 
ad)=x-+Pcost, OR 
Wegen (27) muß 
«—B>1 
gelten; (28) ist erfüllt, wenn man 


ge +2 +24 AT (22 +2+176) 


wählt. Aus diesen Ungleichungen kann man nach 
Vorgabe von « die Größe ß abschätzen. Für « — 2? 
ist z.B. PS 0,21 zulässig. 
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Thermo-Elastic Stresses in an Infinite Slab 
Resting on an Insulating Plane and Having an 
Isolated Heat-Source at the Upper Surface 


= 


Introduction 


The object of this paper is to determine thermo- 
elastic stresses in an infinite slab which is rigidly 
fixed to an insulating surface at the lower base and 
has a point-heat-source at the upper face. Since the 
solution is symmetrical, the method employed by 
SHARMA (1956) has been utilised. 


Method of Solution 


In a steady state of temperature it is known that 
temperature 7’ satisfies the equation. 


a EN en 


Following SHARMA we can deduce the modified 
BELTRAMI equations 


2 
Pot 5a +sETI=0| 
| . - (&)h 
1 2 
er Th 
(3), 


where 
v = PoıIsson’s ratio, 
E = Youn@’s modulus, 
x = Coefficient of linear thermal expansion, 
0=0,+ oy-+ 0, 
and 
PD ,., . an dan 


Solving equations involving Try, Ty; and o, we get 


1 ö 
I = 2(1 9° 8% ee: 
1 ö 
ie een ri hin (5); 
—_— 1 2 
7 Set La N BRRSLE Da ) 


where ®,, Ö,and ®, are harmonic functions to be 
determined by the boundary conditions. 

Let us suppose that the upper surface of the iso- 
tropic elastic slab is given by the plane z = 0, the 
axis of 2 being drawn into the body and the lower 
boundary is the plane 2 = d. 

The temperature is symmetrical about the z-axis 
and hence in steady state it satisfies the equation 


oa I a Rt 


Fe er 
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Also the slab rests on an insulating surface, there- The first two equations give 
fore, when 2=d, B 
sp Tr; = — [m[—mz(Bsinhm2—A cosh m 2) 
Er AN Fr (7). 0 


Equations (6) and (7) are clearly satisfied if we take 


oo 


T= [ M cosh m (e— d) J,(m r) dm (8), 
{) 


where M is an arbitrary function of m and will be 
determined by the boundary conditions. J,(mr) is 
the Bezsser’s function of the first kind and of O-th. 
order. 


Now since the surface of the slab is free we can put 
equations (5) as 
l 


[7) 
ee 


a2 | AEiEh nr nn ln r)lde, 
0 


1 ö 
eg 
> : 0), 
1 j Asinh a [YJ,(m r)] dm , 
0 


1 ö 


aa lteen 


0, = 
= 
+ J B m sinh m z J,(m r) dm 


0 


where A and B are arbitrary functions of m. These 
relations make 775 = Tyz = 05 = 0 ontheplanez = 0. 


Using equations of equilibrium 


we have 

2{+=E7]=2( +») 
oo 

x fr [B cosh m z — A sinh m 2] J,(mr) dm (10), 
0 


which gives d +a@ ET]. 
Substituting in equations (9) we get, 


co 
In= | [—- mz (Bsinhm2— A coshmz) + Asinh mz] 
0 


x = Y,(m r)] dm . 


oo 
BT lan Zoom) Asinhma) (u). 
0 


. 3, Walm r)]dm. 


oo 


%= (} [- m22 (B cosh m2— A sinh m 2) 
Ö 


+ Bm sinh m z] J,(m r) dm 


Ur = — 


+ Asinh mz] Jı(mr) dm. . (12). 
The displacements u, and w are given by 
oo 
j A [2 (1—») cosh mz + mzsinhmz] 
{0 
—B [1 — 2») sinh m z + m z cosh m z] 


nd cosh m (2 — a) au r) dm, 


1-+» 
E 


m 
ES (13). 
el [4 m cosh ms — (1 — 29) sinh ma] 
; 


+ B[2(1—v) cosh mz — mzsinh m z] 
ee 


sinh m (2 — 2) Jo(m r) dm 


To find the stress components or and o, we use the 
relations 


Our U Elle), 
ör a 


Substituting values of 0, o, and u, from the rela- 
tions (10), (11) and (13) we have 


oo 


Ht+og= raw + 2v)sinhmz + mzcosh m2] 
{) 
— A[2(1+v)coshmz-+ mzsinhmz] 
«EM 
— cosh m (<—d) Im r)dm. 
e\ (15), 
= | mia [2(1—v) coshmz + mzsinh m 2] 


Ö 
— B[1— 2»)sinhmz + mzcoshmz] 


ie a cosh m (2 — .) J,(mr) dm 


which give or and o,. 
Now from the boundary condition that the lower 
surface is fixed, we have 


= WI Orewhenzi dar (16): 
From (13) and (16), 
A [2(1—v) cosh m d + mdsinh m d] 
— B[(1—2»)sinh md + mdcosh m d] 
BEN Tun, 


Mm 


A [m d cosh m d — (1 — 2») sinh m d] 

+ B[2 (1—v)coshmd— mdsinhmd] = 0 

The temperature distribution is supposed to be 
such that there is a concentrated heat nucleus at the 
origin so that 


oo 


[2er Te)dr=T, 
0 


or 7, ör) 
1 EM EU ® 
T(r) = mr’ 


where ö(r) is DirAc delta function, 


=, [rum amam IS (18), 
0 


Equation (8) is 
T(r) = [ M cosh mie —d) Ju(m r) dm . 
0 


Comparing we have 


M u 


ee el LO). 
2% 7 cosh md 


Putting v = .25 and 2. — kin equations (17) 


we have r; 
A [3 cosh md + 2m dsinh md] 
— B [sinh md + 2 md cosh md] 
— ksech md = 0,\ (20). 
A [2 md cosh m d — sinh md] : 
+ B[3 cosh m d — 2 md sinh md] = 0 
Solving equations (20) we have, 
k (3 — 2 md tanh md) 
1-+ 4 m?d? + 8 cosh?md ’ 


k (tanh md — 2 md) 
ET + 4 m? d? + 8 cosh? md 


Az 


21), 
B 


A, B being known, allthe stress-components Tyz, Tyz» 
co, are given by equations (11). 
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From equations (11), 
oo 
er a [— m? z (B cosh m z — A sinh m z) 
ö 


+ m Bssinh m 2] J,(m r) dm . 


a 

increase the thickness of the slab, the maxi- 
es of o,/k at the lower base decreases in 
accordance with the figure given below. 


Ve 
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BUCHBESPRECHUNGEN 


Die Umstellung auf das Internationale 
Einheitensystem in Mechanik und Wärme- 
technik. Herausg. H. W. Hahnemann. (Ingenieur- 
wissen 4). 1268. m. 17 Tafeln. Düsseldorf 1959. 
VDI-Verlag. Preis brosch. 9,80 DM. 


Mit eindringlichen Worten und in temperament- 
voller Weise, seine Ausführungen durch eindrucksvolle 
Beispiele belegend, unternimmt der Verfasser mit 
dieser kleinen Schrift einen Feldzug zugunsten der 
allgemeinen Benutzung von ‚Größengleichungen und 
des internationalen Einheitensystems in technisch- 
wissenschaftlichen Veröffentlichungen. Die Vorteile 
der Größengleichungen werden den Nachteilen der 
Zahlenwertgleichungen gegenübergestellt; ein beson- 
deres Kapitel ist dem verwirrenden Durcheinander 
gewidmet, das durch die Wörter ‚Kilogramm und 
Kilopond“ hinreichend gekennzeichnet ist. Das Ver- 
wenden der Masse als Bezugsgröße bei der Aufstellung 
spezifischer Größen wird begründet und sodann das 
internationale Einheitensystem MKSA nach Giorcı 
bzw. das vollständige MKSAKC-System besprochen. 
Beispiele vor allem aus Mechanik und Thermodyna- 
mik (spezifische Wärme, individuelle Gaskonstante, 
Reynoroszahl, Eigenfrequenz usw.) ergänzen und er- 
läutern die Ausführungen. Insbesondere wird der un- 
richtige Gebrauch der Fallbeschleunigung erörtert. 
Die eingefügten Tafeln mit Übersichten über die ver- 
schiedenen Einheiten und ihre Definition und mit 
Umrechnungszahlen sind wertvoll und nützlich. Es 
wäre schön, wenn dieser Schrift der beabsichtigte Er- 
folg beschieden wäre. Jedoch seien dem Referenten 


an dieser Stelle noch einige zusätzliche Bemerkungen 
erlaubt. 


Jeder Ingenieur ist heute daran gewöhnt, mit di- 
mensionslosen Kennzahlen zu arbeiten, ohne daß er 
diesen Zahlen eine besondere Bezeichnung für .die 


Einheit 1 hinzufügt. Warum soll dies ausgerechnet 
beim Bogenmaß eines Winkels geschehen ? Man sollte 
sich bemühen, die Einheit 1 rad wieder auszurotten 
und durch die Zahl 1 zu ersetzen, weil man sonst in 
Schwierigkeiten gerät. Was soll man mit der Einheit 
lrad machen, wenn man beispielsweise die Nähe- 
rungsformel sinx x verwendet? Soll man den 
Sinuswert ebenfalls in rad ausdrücken, weil es mit 
dem Winkel auf der rechten Seite geschieht? Soll 
auch der Einfluß eines in rad gemessenen Winkel- 
fehlers öx auf den Wert von sin x, wenn man ihn nach 
der Formel ösin« = c0os«& - öx berechnet, in rad an- 
gegeben werden? Diese Unzulänglichkeiten kann 
man nur umgehen, wenn man alle Winkeleinheiten als 
Zahlen auffaßt: 1° = n/180, IL = n/2, 17 =r/3200, ... 
Dann ist sin x = « eine (genäherte) Größengleichung, 
bei der es nicht darauf ankommt, in welchem Maß « 
ausgedrückt wird, z. B. 


IT 


=,19-= 


In Tafel 4 der vorliegenden Schrift ist in dieser Hin- 
sicht sowieso eine Inkonsequenz enthalten. Dort ist 
zwar für den Winkel die Einheit rad, aber für die 
Winkelgeschwindigkeit und die Winkelbeschleunigung 
bereits 1/s’ bzw. 1/s? statt rad/s bzw. rad/s? ange- 
geben! 

Man sollte im übrigen überhaupt darauf achten, 
daß hinter den Symbolen der mathematischen Funk- 
tionen wie sin, cos, lg, ... grundsätzlich nur reine 
Zahlen als Argumentwerte vorkommen. Es gibt dann 
z. B. zwar den Logarithmus eines Druckverhältnisses 
In(p/po), aber es gibt nicht den Ausdruck In p — In p,, 
weil die Größen Inp und Inp, einzeln keinen Sinn 
haben. 


Dresden 


sin 10° 10° 0,175. 


H. HEInRıcH 


si daux. 
Gauthier-Villars. 


de Einleitung zu seiner Arbeit umreißt Verf. 
, die Stellung, Bedeutung und Entwicklung des 
ietes der Nichtlinearen Mechanik sowie die beiden 


 topologischen und die quantitativen, die auf Nähe- 
_ rungslösungen hinzielen und in Reihenentwicklungen 
bestehen. Verf. beschränkt sich auf die letzteren und 
behandelt drei Kapitel: Freie Schwingungen nicht- 
linearer Systeme, Sinusförmige Erregung nichtlinearer 
Systeme und Zufallsstörungen in solchen Systemen. 
Er verwendet eine knappe Darstellungsweise, die 

physikalische Gesichtspunkte hervorhebt; einige ma- 
thematische Ergänzungen folgen am Schluß des Hef- 
tes. Im ersten, umfangreichsten Kapitel (55 8.) be- 


aktiven Element (Energiequelle) und einem passiven 
' Element (Übertragungsglied) zusammensetzt, und 
‚erklärt für dieses den komplexen Frequenzgang. Am 
Beispiel des Röhrengenerators zeigt er die Erscheinung 
der konservativen, gedämpften und aufschaukelnden 
Schwingungen; um festzustellen, bei welchem Para- 
meterwert die eine oder andere Erscheinung auftritt, 
benutzt er das graphische Verfahren von Nyauıst, 
das er später auf der Grundlage der harmonischen Li- 
nearisierung für nichtlineare Schwinger ausbaut, um 
- Frequenzen und Amplituden zu berechnen und die 
-  Stabilitätsfrage zu beantworten. Auch noch andere 
Methoden werden erwähnt, und zur Erläuterung wer- 

den die’ einfachsten Gleichungen der Nichtlinearen 
Mechanik (von DUFFING, VAN DER PoL, RAYLEIGH) 
herangezogen. Im 2. Kapitel (25 8.) findet man zu- 
erst die beinahe schon klassischen Untersuchungen 
über das Resonanzdiagramm bei der Durrınsschen 
Gleichung und anschließend interessante Studien mit 
Hilfe der harmonischen Linearisierung. In diesem 

-  $inne wird auch eine Stabilitätstheorie entwickelt. 
- Das 3. Kapitel (35 S.) befaßt sich mit einem neueren 
= Gegenstand, dem Einfluß von Zufallsstörungen. Verf. 
A ‚erörtert erst die Reaktion eines Systems, wenn Ein- 
- 


\ 


Fo. 


heitsstöße als Störungen wirken, führt dann das 
„weiße Rauschen‘‘ auf diese zurück, und bestimmt 
die entsprechenden Amplituden- und Frequenz- 
schwankungen bei den Schwingungen. Alle Überle- 
gungen setzen voraus, daß die Nichtlinearitäten hin- 
_ reichend, schwach sind, so daß die stationären Zu- 
stände der Systeme quasiharmonische Schwingungen 
sind. 
Berlin 


K. Sehütte, Beweistheorie. XI + 355 8. Berlin/ 
Göttingen/Heidelberg 1960. Springer-Verlag. Preis 
geb. 48,— DM. 

Die Hınsertsche Beweistheorie ist entstanden aus 
dem Bestreben, die intuitionsitischen Bedenken gegen 
gewisse Schlußweisen der klassischen Mathematik zu 
entkräften. Diese Bedenken richteten sich vor allem 
gegen eine uneingeschränkte Anwendung des ‚‚ter- 
tium non datur‘‘, vermöge dessen in der klassischen 
Mathematik zum Beispiel eine Existenzaussage be- 
reits dann als bewiesen gilt, wenn die Annahme der 
Nichtexistenz zum Widerspruch geführt ist (ein 
Intuitionist sieht eine Existenzaussage nur dann als 
bewiesen an, wenn ein als existent behauptetes Ob- 
jekt wirklich konstruiert wurde). Hierzu gab HILBERT 
der klassischen Mathematik eine neue Deutung, indem 
er die jeweils betrachtete mathematische Theorie einer 
'Formalisierung unterwarf und die Beweiszusammen- 
hänge der klassischen Mathematik als rein formale 
Zusammenhänge zwischen den die mathematischen 
Sätze symbolisierenden Formeln deutete — ein 
mathematischer Beweis wurde also zu einer nach be- 
stimmten formalen Gesetzen gebildeten endlichen 
Folge von Formelfiguren, wobei diese Gesetze in 


R. Reıszıe 


ee en = 


= 


Les 


klassischen Mathematik (unter Einschluß also z. B. 


wichtigen Untersuchungsmethoden: die qualitativen 


 trachtet Verf. einen Schwinger, der sich aus einem 


A % Bi, 
> 


Ki 
des tertium non datur) wiedergeben. Bei dieser Auf- 


fassung äußert sich die Frage nach der Sinnhaftigkeit 
der klassischen Schlußweisen im Problem der formalen 


Widerspruchsfreiheit des betrachteten formalen Sy- 
stems, d. h. ob man nicht etwa mit den formalen Um- 
formungsregeln aus den betrachteten Axiomen eine 
gewisse Formel und ihr Negat erhalten kann. Die 


Aufgabe der Beweistheorie besteht nun darin, diesen 


Widerspruchsfreiheitsbeweis mit ‚„‚finiten‘‘ Mitteln zu 
führen, also mit solchen, wie sie von den Intuitionisten 
anerkannt werden. Hierzu reichen nun zwar nach 
einem berühmten Resultat von GÖDEL die engsten 
finiten Mittel nicht aus, jedoch gelingt z. B. der Wider- 
spruchsfreiheitsbeweis für die elementare Arithmetik 
der natürlichen Zahlen (nach AcKERMANN und 
GENTZEN) mit transfiniten Hilfsmitteln, die in einem 
ganz bestimmten Sinne durchaus noch als konstruktiv 
anzusehen sind (transfinite Induktion bis zur ersten 
e-Zahl). 

Der vorliegende Band gibt nun in bemerkenswerter 
Vollständigkeit die Resultate der Untersuchungen zur 
Beweistheorie, die in den mehr als zwanzig Jahren er- 
zielt wurden, die seit dem Erscheinen der ‚„‚Grund- 
lagen der Mathematik‘ von HıLBERT und BERNAYS 
vergangen sind. 

‚In einem ersten Teil werden zunächst einerseits die 
Hilfsmittel aus der Logik und andererseits die Hilfs- 
mittel aus der Theorie der (konstruktiven) Ordinal- 
zahlen entwickelt. Die Logik wird dabei in einer dem 
GENTzEnschen Sequenzenschließen verwandten Form 
dargestellt, die auch für sich von gewissem Interesse 
ist. Bezüglich dieses ersten Teiles ist vor allem zu be- 
merken, daß der Autor bewußt darauf verzichtet hat, 
eine umfassende Darstellung der Logik oder der 
Ordinalzahltheorie zu geben, sondern sich auf die für 
die späteren Untersuchungen notwendigen Dinge be- 
schränkt — ohne daß dadurch die für eine Mono- 
graphie notwendige Systematik verloren geht. 

In einem zweiten Teil werden dann zunächst die 
klassischen Resultate über die beweistheoretischen 
Fragen der elementaren Zahlentheorie behandelt. Im 
Anschluß daran wird ein finiter Widerspruchsfreiheits- 
beweis für das System der elementaren Zahlentheorie 
mit unendlicher Induktion gegeben. Den Abschluß 
bilden Untersuchungen über die in der elementaren 
Zahlentheorie herleitbaren Formen der transfiniten 
Induktion. 

In einem dritten Teil werden schließlich beweis- 
theoretische Resultate zu verschiedenen formalen 
Systemen der Mengenlehre und der Analysis ent- 
wickelt, und damit wird der Anschluß an die im gegen- 
wärtigen Zeitpunkt aktuellen Probleme errreicht. 

Die vorliegende Monographie kann jedem Mathe- 
matiker, der sich mit den Fragen der Beweistheorie 
vertraut machen möchte, wärmstens empfohlen 
werden. Als Lehrbuch kommt es naturgemäß allen- 
falls für Studenten der letzten Semester in Frage. 

Greifswald G. ASsER 


Prof. Dr. S. Valentiner, Vektoren und Matrizen‘ 
(Sammlung Göschen, Band 354/354a.) 2. Aufl. 202 8. 
Berlin 1960. Walter de Gruyter & Co. Preis brosch. 
5,80 DM. 

Das aus der ‚„‚Vektoranalysis‘‘ des nämlichen Ver- 
fassers hervorgegangene Bändchen unterscheidet sich 
von der voraufgehenden Auflage gleichen Titels nur in 
ganz wenigen Punkten. Die $11—13, in denen 
Vektor- und Skalarprodukte von Faktoren, die selbst 
schon Vektorprodukte sind, behandelt werden, haben 
einige geringfügige mathematische Änderungen er- 
fahren; im $14 wird der Frage der Division von 
Vektoren auf einer halben Seite Aufmerksamkeit ge- 
schenkt. Im übrigen vgl. die Besprechung der letzten 


Auflage (ZAMM 38 (1958), S. 487—488). 
Dresden L. BITTNER 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 
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Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten, 


tatist 


rn 


iewin Automatic Programming. R. Kurth, Axiomatics of Classical 8 
ee A., B.Se. x +3208. cal Mechanics. X + 180 8. Oxford/Lon 
Oxford / London /New York /Paris 1960. Pergamon York/Paris 1960. Pergamon Press, 
Press. Preis geb. 63 s. net. 45 s. net. 


J.Haag et R.Chaleat, Problömes de theorie H. Parkus, Mechanik der festen er. VOL 
gönerale des oscillations et de chronomötric. 1264 S. m. 191 Abb. Wien 1960. Spri Verlag. x 
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Dr 


412 S. m. 66 Abb. Paris 1960. Gauthier-Villars. Preis geb. DM 29,—. 


Preis brosch. $ 8,00. 


U. Bornitz und A. Schubert, Siebenstellige Log- 
arithmen der Fakultäten der Zahlen 0 bis 
12000. 518. Leipzig 1960. B. G. Teubner Verlagsge- 
sellschaft.. Preis geb. DM 10,30. 


I. I. Eterman, Analogue Computers. IX -+264S8. 
m. 123 Abb. u. 38 Tafeln. Oxford/London/New York/ 
Paris 1960. Pergamon Press. Preis geb. 50 s net. 


B. Spain (Head of the Mathematics Department, 
Sir John Cass College), Analytical Quadrics. 
IX + 1358. m. 23 Abb. Oxford/London/New York/ 
Paris 1960. Pergamon Press. Preis geb. 30 s net. 


H. Wolf (Direktor des Instituts für theoretische 
Geodäsie, Bonn), Ausgleichsrechnung nach der 
Methode der kleinsten Quadrate. 488. Ham- 
burg 1960. Hanseatische Verlagsanstalt GmbH. Preis 
brosch. DM 6,80. 


H. E. Salzer, D. €. Shoultz, and E.P. Thompson, 
Tables of Osculatory Integration Üoeffi- 
cients. 438. San Diego 1960. Convair (Astronau- 
tics), Division of General Dynamics Corporation. 


R. Schatten, Norm Ideals of Completely Con- 
tinuous Operators. (Ergebnisse der Mathematik 
und ihrer Grenzgebiete, Heft 27). VIII+81S. 
Berlin/Göttingen/Heidelberg 1960. Springer-Verlag. 
Preis brosch. DM 23,60. 


I. Babuska — K. Rektorys — F. Vyeichlo, Mathe- 
matische Elastizitätstheorie der ebenen Pro- 
bleme. VIII - 478 S. m. zahlr. Abb. Berlin 1960. 
Akademie-Verlag. Preis geb. DM 96,—. 


D. Laugwitz, Differentialgeometrie. 183 S. m. 
44 Abb. Stuttgart 1960. B.G. Teubner Verlags- 
gesellschaft. Preis geb. DM 24,60. 


A.D. Smirnov, Tables of Airy Functions and 
SpecialConfluent Hypergeometrice Funetions. 
(Mathematical Tables Series, Volume 7). VII + 260 8. 
Oxford/London/New York/Paris 1960. Pergamon 
Press. Preis geb. £5 net. 


W. Kämmerer, Ziffernrechenautomaten. 2. un- 
veränderte Auflage. VIII -+ 3038. m. 156 Abb. 
Berlin 1960. Akademie-Verlag. Preis geb. DM 29,—. 


M. R. Redwood (Lecturer in Rlectrical Engineering, 
Quenn Mary College, Univ. of London), Mechanical 
Waveguides. The propagation of acoustie and ultra- 
sonic waves in fluid and solids with boundaries. VIII 
+ 300 8. Oxford/London/New York/Paris 1960. Perga- 
mon Press. Preis geb. 50 s. net. 


Archive for History of Exact Sciences. Vol.l, 
Number 1. Edited by C. Truesdell. 106 8. Berlin/ 
Göttingen/Heidelberg 1960. Springer-Verlag. Preis 
brosch. DM 19,60. 


J. Tsehauner, Einführung in die Theorie der 
Abtastsysteme. 185 S. m. 81 Abb. u. 5 Oszillo- 
grammen. München 1960. Verlag R. Oldenbourg. 
Preis brosch. DM 32,—. 


6. Otto, Moderne Leichtbautechnik, Band]. 
Rumpf, Vollwand- und Schalenbauweise. XII +284 S. 
m. 201 Abb. Braunschweig 1960. Friedr. Vieweg 
& Sohn. Preis brosch. DM 24,80. 


W. Häussler, Das Mollier-i, <-Diagramm für 
feuchte Luft und seine technischen Anwen- 


dungen. XII + 167 S. m. 92 Abb. u. 9 Beilage- _ 


tafeln. Dresden und Leipzig 1960. Verlag Theodor 
Steinkopff. Preis geb. DM 18,30. 


R. M. L. Baker, Jr., and M. W. Makemson (Depart- 
ment ‘of Astronomy, University of California), An 
Introductionto Astrodynamies. XIV + 3588. 
m. 47 Abb. New York and London 1960. Academic 
Press. Preis geb. $ 7,50. 


E. B. Dynkin, Theory of Markov Processes. 
Translated from the Russian by D. D. Brown, M. A. 
IX + 210 S. Oxford/London/New York/Paris 1960. 
Pergamon Press Ltd. Preis geb. 60 s. net. 


Russian Mathematical Surveys. Vol. 15, 
Numb. 1 (Uspekhbi Matematicheskikh Nauk). 130 S. 
London 1960. Cleaver-Hume Press Ltd. Preis brosch. 
single number £ 2. 


BERICHTIGUNGEN 


Zu J. Brilla, Einige gemischte Randbedingungen 
anisotroper Platten, ZAMM 40 (1960), S. 472—482: 


Nach Mitteilung des Verfassers lautet der Nenner 
auf der rechten Seite der Formel (5.11) 2 ri und der 
Nenner auf der rechten Seite von Formel (5.19) 4 x A. 


Zu W. Franz, Zur Asymptotik der Zylinderfunk- 
tionen und ihrer Nullstellen, ZAMM 40 (1960), 8. 392: 


In Formel (12 b) ist ein Druckfehler unterlaufen. 
Im ersten Glied der dritten Zeile lautet der erste 
Faktor im Nenner 125 (anstatt 152). 
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SIEGFRIED KÄSTNER 
- "Vektoren, Tensoren, Spinoren 
Eine Einführung in den Tensorkalkül unter Berücksichtigung der physikalischen Anwendung 


1960. XII, 308 Seiten — 27 Abbildungen — gr. 8° — Ganzleinen DM 38,— 


Will man tiefer in die Gesetze der Physik eindringen, so muß man sich mit dem Wesen und 
den Eigenschaften von Tensoren vertraut machen, da sämtliche Größen der Physik Tensoren 


sind. 


Das entsprechende mathematische Werkzeug ist der allgemeine Tensorkalkül, der unter Ver- 
wendung einer der Natur der Tensoren angepaßten Schreibweise eine übersichtliche Formu- 


lierung sämtlicher Tensorbeziehungen gestattet. 


Der Verfasser geht vom geometrisch anschaulichen Vektorkalkül im dreidimensionalen 
euklidischen Raum aus und macht den Leser schrittweise mit den Grundlagen des allgemeinen 
Tensorkalküls bekannt. Dabei werden die eingeführten Begriffe an Hand zahlreicher Rechen- 


beispiele und physikalischer Beispiele erläutert. 


Schließlich werden auch Tensoren im N-dimensionalen euklidischen bzw. pseudo-euklidischen 
Raum und im N-dimensionalen RIEMANNschen Raum behandelt, die im Rahmen der spe- 
ziellen und allgemeinen Relativitätstheorie sowie der Quantentheorie eine maßgebende Rolle 


spielen. 


Bestellungen durch eine Buchhandlung erbeten 
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Vorschläge zur \ Vi 
in die Berliner A Ak 


(Abhandlungen der Deutschen Akademie der \ 
Physik und Technik, Jahrgang 1960, Heft BR 


1960. 74 Seiten — 4° — DM 7,50 


‚Fachurteile über ieh zur Zeit der I 
in Berufsverhandhungen der Universitäten, 
brieflichen Äußerungen hauptsächtlich in den Wah 
‚Dirichlet war der erste, der in solchen Vorschlägen Pan 
lichkeit und Leistung der zur Wahl in die Berliner / 
geben hat. Nach ihm haben u.a. vor allem ER TE 
reiche und stilistisch vollendete „‚Landationes“ eingereicht. Der We e 
derartiger Vorschläge zur Wahli in die Hoikaer Akademie, intern 

. des vorigen Jahrhunderts, für en a er vor und e 


mit Te Ne über ER NE und ER, der W 
bindet. , 


ERICH KÄHLER 


Innerer und äußerer Differentialkalkül 


(Abhandlungen der Deutschen Akademie der Wissenschaften zu Hecin, Klasse fü 
Physik und Technik, Jahrgang 1960, Heft 4) 


1960. 32 Seiten — 4° — DM 5,50 
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